file  COPY  ADA112093 


VSC-TR-8 1-29 


VFM  DISCRIMINATION 
RESULTS  FROM  A 
TEN  STATION  NETWORK 


J.  M.  Savino 
C.  B.  Archambeau 
J.  F.  Masso 


a  20  c 


SYSTEMS,  SCIENCE  AND  SOFTWARE 

P.  O.  Box  1620 

La  Jolla,  California  92038 


July  1980 


TECHNICAL  REPORT 


APPROVED  FOR  PUBLIC  RELEASE, 
UNLIMITED  DISTRIBUTION 


Copy  available  to  DT1C  does  Uvi 

permit  fully  legible  reproduction 


Monitored  by: 

VELA  Seismological  Center 
312  Montgomery  Street 
Alexandria,  Virginia  22314 


P' 


J*. 


MAR  1  7  1982 


s 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


i 


_ UNCLASSIFIED _ 

SECURITY  CL  ASSlF  iCATIQm  OF  TnlS  PAGE  -»~K»n  Oaia  Enund)  _ __ 

r - ^....rviTiTiAil  Circ  '  READ  INSTRUCTIONS 

REPORT  DOCUMENTATION  PAGE  !  before  completing  form 

I  REPOST  numbs" - 2  SOVT  ACCESSION  NO.!  1  REClPI  En  t-S  C»T»l:5  NUMBER 

VSC-TR-81-29  PVbAllO  £4$  _ 


1  REPOST  NUMBER 


4.  TITLE  'end  Subtitle) 

VFM  DISCRIMINATION  RESULTS  FROM  A 
TEN  STATION  NETWORK 

7  AUTHORED 

J.  M.  Savino 
C.  B.  Archambeau 

J.  F.  Masso  _ _ _ 

9.  PERFORMING  ORGANIZATION  NAME  and  AOORESS 

Systems,  Science  and  Software 

P.  0.  Box  1620 

La  Jolla,  California  92038 


I  5.  type  op  report  4  period  ddvered 

j  Technical  Report 
10/1/79  -  9/30/80 

rr  PERFORMING  3RG.  REPOR"  n.mBER 

|  SSS-R- 80-4566 

I  9.  CONTRACT  OR  GRANT  NUMSEA’ii 


F08606-79-C-0008 


|  ’0.  program  element  projEC'  ~  as* 
AREA  4  RCRK  jnit  NUMBERS 

Program  Code  No.  6H189 
ARPA  Order  No.  2551 


71  CONTROLLING  OFFICE  NAME  ANO  AOORESS  I  12-  REPORT  DATE 

VELA  Seismological  Center  I  July  9  0 _ 

312  Montgomery  Street  j  ”  NCM8ER  °r  "*ses 

Alexandria,  Virginia  22314 _ 1  ,.T: - — 

;  4  MONITORING  AGENCY  NAME  4  AOORESS,//  lllttnri  from  Corure,i,n«  Otlic.i pS.  SECUPI  .  V  .LASS,  ol  thtw  topor 

Unclassified 


Ts«  2ECL  ASSl  FlC  A’lQN  DOWNGRACiNG 
SCHEDULE 


hr  DISTRIBUTION  STATEMENT  rol  this  Report ) 


Approved  for  Public  Release,  Distribution  Unlimited. 


DISTRIBUTION  STATEMENT  'ol  the  ebetrmct  entered  in  3lock  20.  tl  different  from  Report) 


18.  supplementary  notes 


19.  KEY  *OR9S  ' Continue  on  reverse  side  il  necesssrv  end  identify  8v  block  number) 

Seismic  Discrimination 
Narrow  Band  Filtering 
Earth  Noise 

Anelastic  Attenuation  .  .  r  /  .... 


n  /  s 


20.  ABSTRACT  ( Continue  on  reverse  side  If  necessmrv  end  identify  bv  block  number ) 

In  this  report  we  describe  the  results  of  a  discrimination 
experiment  performed  on  a  large  number  of  Eurasian  events.  Our 
objective  in  this  experiments  was  to  identify  each  of  these  events 
as  either  an  earthquake  or  an  underground  explosion  based  on  an 
analysis  of  short-period  P  waves  recorded  at  a  global  network  of 
seismograph  stations.  rIn  all,  there  were>ten  seismograph  station 
involved  in  this  experiment.  -  (continued) 


_ IINPLASSIFIED _ _ 

SECURITY  CLASSIFICATION  OF  TwiS  PAGE  'l*h»n  Daia  I 


DD  ,  1473  EOlTION  OF  ■  NOV  S9  IS  OBSOLETE 


3-jIo  U  '") 

_ UNCLASSIFIED  _ 

SECURITY  CLRSSiriCRTIQH  OF  THIS  ■  ■” — “ — 

20.  ABSTRACT  (continued) 

The>dis crimination  results  obtained  from  this  experiment 

support  the  theoretical  predictions  of  Archambeau,  e|  al.,  1964 
suppui-u  i-i  .  .  fhp  variable  frequency  magnitude  (VFM) 

rto^noSvid^  verfef^cuve  discriminLion  between  Eurasian 
XftSqua/erinfe^oLone  recorded  at  both  teleeeramrc  and  re- 

?  Hicfancps  One  of  the- principal  findings  of  this  stuay 

is  that  the  effectiveness  of  mb(f)  type  discrimination  at  each 
of  the  stations  included  in  this  experiment  is  largely  controlle 

j.sfis.?  sssssJyW: 

earthquakes  and  explosions  than  do  mb(f)  estimates  from 
tion  located  over  low  Q  upper  mantle. 


_ unclassified _ 

SECURITY  classification  or  this  RRGEfWtim  Oar*  Enr .red) 


AFTAC  Project  Authorization  No.  VT/0712/B/PMP 
ARPA  Order  No.  2551,  Program  Code  No.  6H189 
Effective  Date  of  Contract  November  17, 1978 
Contract  Expiration  Date:  November  15, 1981 
Amount  of  Contract  $1,816,437 
Contract  No.  F08606-79-00008 
Principal  Investigator  and  Phone  No. 

Dr.  J.  Theodore  Cherry,  (714)  453-0060 

Project  Scientist  and  Phone  No. 

Mr.  Brian  W.  Barker,  (202)  325-7581 


This  research  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  and  was  monitored  by 
AFTAC/VSC,  Patrick  Air  Force  Base,  Florida  32925,  under  Contract  No. 
F08606-79-C0008. 

The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  Advanced 
Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  U.S.  Government 


_Accesslon  For 

NTIS  GRA&I 
D71C  TAB  ^ 

Un.'-.nnounord 

Ju:  t  i  r  i  c  !,  t  -?r, 


Py .. 

r  i  r- -  • 


Aval.':  i 


Codo3 


Diet 


'ur 


a 


I  Js 

Up 


DD 


TABLE  OF  CONTENTS 


Section 

I. 

II. 


III. 


IV. 


V. 

APPENDIX  A 


Page 


INTRODUCTION .  1 

AUTOMATED  DETECTION  AND  DISCRIMINATION.  ...  5 

2.1  VFM  DISCRIMINATION  PROCEDURE  .  5 

2.1.1  Narrow-Band  Filtering  .  5 

2.1.2  The  Automated  Detection 

Algorithm . 20 

2.1.3  The  Discrimination  Algorithm.  .  .  32 

VFM  DISCRIMINATION  RESULTS  FROM  A  TEN 

STATION  NETWORK  .  47 

3.1  VFM  RESULTS  FROM  TELESEISMIC  STATIONS.  .  47 

3.2  VFM  RESULTS  FROM  STATIONS  AT  REGIONAL 

AND  TELESEISMIC  DISTANCES . 57 

DISCUSSION  AND  CONCLUSIONS . 71 

4.1  INFLUENCE  OF  PATH  AND  SOURCE  EFFECTS 

ON  VFM  DISCRIMINATION . 71 

4.2  CONCLUSIONS . 77 

REFERENCES . 7  8 

:  TABULATION  OF  DISCRIMINATION  PARAMETERS  ...  80 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


3. 


4. 


5. 

6. 

7. 

3. 

9. 

10a. 

10b. 


LIST  OF  ILLUSTRATIONS 

Page 


MARS  flowchart .  6 

Flowchart  indicating  principal  mathematical 
operations  in  the  MARS  program .  9 


Examples  of  the  outputs  of  a  suite  of  4.0 
to  0.5  Hz  narrow  band  filters  applied  to  a 
seismogram  from  a  Kuril  earthquake  recorded 
at  a  station  at  teleseismic  distances . 18 

Examples  of  the  outputs  of  a  suite  of  narrow 


band  filters  applied  to  a  seismogram  of  a 

presumed  explosion  recorded  at  the  same 

station  as  in  Figure  2a . 19 

Typical  tg-f  plane  representation  of  a  time 
series  segment  (0  to  36  seconds)  when  the 
signal- to-noise  ratio  is  low . 22 

A  seismogram  written  at  the  Seismic  Research 
Observatory  (SRO)  in  Chiang  Mai,  Thailand 
(CHTO)  from  a  presumed  explosion  in  eastern 
Kazakhstan . 24 


Seismogram  written  at  Tatalina,  Alaska,  for 

the  same  presumed  explosion  as  in  Figure  4  .  .  25 


Seismogram  recorded  at  a  station  in  Bluff, 

Alaska,  from  an  event  in  the  Kurils . 26 

Seismogram  recorded  at  NORSAR  for  the  same 
Kuril  earthquake  shown  in  Figure  6 . 27 

Example  of  timing  of  multiple  phase  arrivals 
by  MARS . 29 

Spectral  magnitudes,  mfc(f),  computed  at  0.45 
Hz  and  2. 25  Hz . 33 

Spectral  magnitude  estimates  at  fc  *  0.6  Hz 
and  fc  =  5.0  Hz  for  an  event  population  re¬ 
corded  at  the  Oyer  array  in  Norway . 35 

Maximum  filter  amplitudes  with  noise  correc¬ 
tion  for  the  same  event  population  plotted 
in  Figure  10a . 36 


ii 


SYSTEMS.  SCIENCE  AND  SO FTWaH 


LIST  OF  ILLUSTRATIONS  (continued) 


Figure  Page 


11a.  Typical  event  distributions  in  the  ntb(f) 

plane  for  event  data  that  is  not  corrected 

for  noise  contamination . 45 

lib.  Typical  event  distributions  in  the  mfc(f) 

plane  for  noise  corrected  event  data . 46 

12.  VFM  results  for  the  station  located  at  Red 

Lake,  Ontario  (RKON) . 48 

13.  VFM  results  for  the  station  located  at 

Charters  Towers,  Australia  (CTAO) . 51 

14.  VFM  results  for  the  station  located  at 

Albuquerque,  New  Mexico  (ANMO)  .  52 

15a.  VFM  results  for  all  events  recorded  at  the 

station  in  Taipei,  Taiwan  (TATO)  .  54 

15b.  VFM  results  for  the  station  located  at 

Chiang  Mai,  Thailand  (CHTO) . 55 

15c.  VFM  results  for  the  station  located  at 

Zongo  Valley,  Bolivia  (ZOBO)  .  56 

16a.  VFM  results  for  the  station  located  at 

Houlton ,  Maine  (HNME) . 58 

16b.  VFM  results  for  the  center  subarray  at 

LASA  (LAO) . 59 

17.  Map  of  Eurasian  earthquakes  and  presumed 

explosions  located  between  15  and  23 
degrees  from  the  SRO  station  at  Kabul, 

Afghanistan  (KAAO)  .  60 

18a.  VFM  discrimination  results  for  Kabul, 

Afghanistan  (KAAO)  .  61 

18b.  VFM  results  for  teleseismic  events  recorded 

at  KAAO . 63 

19a.  VFM  results  for  regional  events  recorded  at 

ILPA . 64 

19b.  VFM  results  for  teleseismic  events  recorded 

at  ILPA . 65 


in 


SYSTEMS.  SCIENCE  AN O 


LIST  OF  ILLUSTRATIONS  (continued) 


Figure  Page 


20a.  VFM  results  for  the  station  at  Mashhad, 

Iran  (MAIO)  for  events  in  the  distance 

range  18. 7°  <  A  <  28. 7° . 67 

20b.  VFM  results  for  teleseismic  events  re¬ 
corded  at  MAIO . 68 

21.  Teleseismic  P-wave  residuals  based  on  data 

from  the  bulletins  of  the  ISC  for  seismo¬ 
graph  stations  in  the  region  of,  and 
including,  MAIO,  KAAO  and  ILPA . 69 


22a.  VFM  discrimination  results  for  Kabul, 

Afghanistan  (KAAO)  for  the  earthquakes 
and  presumed  explosions  mapped  in  Figure  17. 


These  events  are  located  between  15°  and 

23°  from  KAAO . 74 

22b.  VFM  results  for  the  station  located  at  Red 

Lake,  Ontario  (RKON) . 75 


iv 


SYSTEMS .  SCIENCE  ANO  SOFTWARE 


> 


I.  INTRODUCTION 

In  this  report  we  describe  the  results  of  a  discrimina¬ 
tion  experiment  performed  on  a  large  number  of  Eurasian  events. 
Our  objective  in  this  experiment  was  to  identify  each  of  these 
events  as  either  an  earthquake  or  an  underground  explosion 
based  on  an  analysis  of  short-period  P  waves  recorded  at  a 
global  network  of  seismograph  stations.  In  all,  there  were 
ten  seismograph  stations  involved  in  this  experiment.  Results 
from  related  studies  are  presented  in  Savino,  et  al.  (1979) 
and  Savino,  et  al.  (1980) . 

The  discrimination  technique  that  was  used  in  this 
experiment  is  referred  to  as  the  variable  frequency  magnitude 
(VFM)  technique.  A  fairly  detailed  formulation  of  the  VFM 
approach  and  several  examples  of  its  application  as  an  auto¬ 
mated  detection-discrimination  technique  will  be  discussed  in 
later  sections  of  this  report.  Briefly,  the  central  feature 
of  the  VFM  method  is  the  use  of  a  set,  or  "comb,"  of  Gaussian 
shaped  narrow  band  filters  to  decompose  a  digitial  time  series, 
consisting  of  signal  plus  noise,  into  a  set  of  quasi-harmonic 
modulated  signals.  The  Gaussian  form  of  the  filters  ensures 
optimal  time  and  frequency  resolution  within  the  constraints 
imposed  by  the  sampling  theorem  or  uncertainty  principle  (by 
contrast,  purely  harmonic  components,  or  Fourier  spectra,  con¬ 
tain  only  frequency  information) .  The  filter  center  frequencies 
are  selected  to  span  the  frequency  range  of  the  expected  signal 
spectra.  For  each  filter  a  time  domain  envelope  function  is 
constructed  from  the  filtered  and  quadrature  signals.  The  time 
of  the  maximum  of  a  particular  envelope  function  is  the  group 
or  energy  arrival  time,  tg,  of  "signal"  with  a  frequency  near 
the  filter  center  frequency.  The  amplitude  cf  the  envelope 
function  at  time,  tg,  is  the  spectral  amplitude,  Ag(f),  of  the 
signal  energy  arrival  at  the  filter  center  frequency. 
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For  identification  of  a  signal  as  to  the  originating 

event  type,  the  procedure  is  simply  to  compute  the  magnitude 

m.  (f)  («  log  A  (f))  and  construct  m.  (f)  planes  wherein  event 
o  g  d 

m^(f)  values  at  a  high  frequency  (e.g.,  >  2.0  Hz)  are  plotted 
versus  m^f)  values  at  a  lower  frequency  (e.g.,  0.5  Hz). 
Observed  signals  can  then  be  characterized  by  the  location  of 
their  m^(f)  values  in  these  planes.  It  is  predicted  theoreti¬ 
cally  (Archambeau,  et  al . ,  1974)  that  ordinary  earthquake  sig¬ 
nals  fall  in  one  particular  area  within  such  planes  and  that 
explosions  (of  all  types)  fall  within  a  distinctly  separate 
area. 


The  m^f)  results  obtained  from  this  experiment  support 
the  theoretical  predictions  of  Archambeau,  et  al_.  (1974)  and 
demonstrate  that  the  VFM  approach  provides  a  very  effective 
discriminant  between  Eurasian  earthquakes  and  explosions  re¬ 
corded  at  both  teleseismic  and  regional  distances.  One  of 
the  principal  findings  of  this  study,  however,  is  that  the 
effectiveness  of  m^(f)  type  discrimination  at  each  of  the 
stations  included  in  this  experiment  is  largely  controlled  by 
the  attenuation  properties  of  the  upper  mantle  beneath  a  par¬ 
ticular  station.  For  instance,  m^Cf)  estimates  from  data 
recorded  at  a  station  overlying  high  Q  upper  mantle  (as 
evidenced  either  from  the  location  of  the  station  in  a  shield 
region  and/or  a  region  characterized  by  large  negative  P-wave 
travel-time  residuals)  consistently  yield  a  greater  degree  of 
separation  of  earthquakes  and  explosions  than  do  m^ (f )  esti¬ 
mates  from  a  station  located  over  low  Q  upper  mantle.  Based 
on  this  result,  we  devised  a  fairly  simple  station  classifica¬ 
tion  scheme  (Savino,  et  al,  1979  and  Savino,  at  al. ,  1980). 
Accordingly,  stations  were  specified  as  being  a  Class  I,  II 
or  III  type  depending  upon  discrimination  performance  and 
evidence  of  upper  mantle  attenuation  properties  —  a  Class  I 
station  providing  the  greatest  separation  of  earthquake  and 
explosion  populations,  a  Class  III  station  the  least.  In 
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terms  of  this  classification  scheme,  the  principal  conclusions 
resulting  from  this  experiment  are  the  following: 

•  VFM  discrimination  can  be  achieved  to  im  ^  4 
at  Class  I  regional  stations. 

•  VFM  discrimination  based  on  both  regional  and 
teleseismic  stations  is  approximately  95  per¬ 
cent  effective  for  all  the  events  in  this 
experiment. 

•  The  VFM  technique  is  an  effective  discrimi¬ 
nant  principally  because  of  the  enriched 
high  frequency  source  spectra  of  explosions 
relative  to  earthquakes  with  the  same  spec¬ 
tral  level  near  two  seconds  period. 

The  VFM  discrimination  procedure  is  embodied  in  the 
MARS  (Multiple  Arrival  Recognition  System)  computer  program. 
During  the  course  of  this  discrimination  experiment,  and  under 
support  from  both  the  Vela  Seismological  Center  (VSC)  and  the 
U.  S.  Arms  Control  and  Disarmament  Agency  (ACDA) ,  an  automated 
detection-discrimination  algorithm  based  on  the  narrow  band 
filtering  technique  was  developed  and  added  to  MARS.  The 
approach  adopted  in  the  formulation  of  the  detection  algorithm 
(Masso,  et  al. ,  1979)  is  one  of  pattern  recognition.  That  is, 
based  on  the  output  from  a  comb  of  narrow  band  filters  applied 
to  a  digital  seismogram  the  algorithm  searches  for  patterns  in 
the  signal  and  noise  information  as  it  is  expressed  in  the 
t  -f  plane  (the  group  arrival  time  versus  frequency  plane) . 

In  the  case  of  a  body  wave,  the  pattern  to  be  searched  for  is 
that  of  a  (nearly)  undispersed  signal  with  spectral  amplitudes 
significantly  above  background  noise.  Thus,  in  the  t  -f  plane 
the  algorithm  searches  for  a  straight  line-up  in  time  of 
envelope  maxima  over  a  band  of  frequencies.  Included  in  the 
algorithm  is  a  somewhat  unique  approach  to  the  treatment  of 
"noise"  contamination  of  detected  signals.  All  the  event 
seismograms  included  in  the  discrimination  experiment  were 
processed  with  this  automated  version  of  MARS. 
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The  organization  of  this  report  is  as  follows.  In 
Section  II  we  describe  the  narrow  band  filtering  technique 
and  the  automated  detection-discrimination  algorithms  that 
comprise  the  VFM  technique  used  in  this  experiment. 

In  Section  III  we  present  the  VFM  discrimination  results 
for  each  of  the  ten  stations  included  in  this  phase  of  the 
experiment.  A  tabulation  of  all  the  VFM  discrimination  param¬ 
eters  is  given  in  Appendix  A  for  each  station.  Multistation 
estimates  of  event  types  (explosion  or  earthquake)  based  on 
the  ten  stations  included  in  this  phase  of  the  experiment  and 
a  number  of  other  stations  are  discussed  in  Savino,  et  al. 
(1980). 
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II.  AUTOMATED  DETECTION  AND  DISCRIMINATION 


2.1  VFM  DISCRIMINATION  PROCEDURE 

The  variable  frequency  magnitude  (VFM)  discriminant,  as 
originally  proposed  by  Archambeau,  et  al.  (1974)  and  subse¬ 
quently  developed  and  applied  by  Savino  and  Archambeau  (1974) , 
is  the  basic  technique  that  was  used  to  identify  the  Eurasian 
events  included  in  this  experiment.  While  the  VFM  approach  is 
completely  general  in  its  application,  that  is,  it  can  be 
applied  to  multicomponent  long  and  short  period  data,  the 
emphasis  here  is  on  short  period  vertical-component  P-wave 
seismograms. 

2.1.1  Narrow  Band  Filtering 

The  central  feature  of  the  VFM  signal  analysis  approach 
is  the  use  of  narrow  band  frequency  filters  to  break  up  or 
decompose  a  time  series  consisting  of  signal-plus-noise  into 
a  set  of  quasi-harmonic  modulated  "signals."  This  set  of 
filtered  signals,  one  for  each  filter  center  frequency,  can 
then  be  used  to  determine  the  energy  arrival  time  (the  group 
time,  tg)  and  amplitude  of  the  signal  for  each  center  frequency 
by  analysis  of  the  time  modulation  of  the  filter  outputs. 

Over  the  past  several  years,  the  VFM  signal  analysis 
procedure  has  undergone  several  stages  of  refinement.  The 
automated  procedure  that  was  used  during  this  experiment  is 
embodied  in  a  large  versatile  computer  program  referred  to  as 
the  Multiple  Arrival  Recognition  System  (MARS) .  Figure  1 
summarizes  the  flow  of  operations  in  the  MARS  program:  Figure 
la  provides  a  verbal  outline;  Figure  lb  presents  the  key 
mathematical  operations  performed  in  this  program.  Although 
not  exercised  on  the  data  set  in  this  experiment,  the  polariza¬ 
tion  filtering  techniques  that  have  been  implemented  in  MARS 
are  included. 
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LIST: 

1.  CONTROL,  PLOT,  & 
REPORT  OPTIONS 

2.  INPUT  DATA 

3.  NBF  PARAMETERS 


READ  DESIRED  EVENT, 
SINGLE  COMPONENT  TIME  SERIES 
LOOP  SIGNAL  AND  NOISE 

OVER  _ _ _ 

ADDITIONAL  } 

COMPONENTS  ? 


PLOT  INSTRUMENT 
TRANSFER  FUNCTION 


SCALE,  DEMEAN,  DETREND 
AND  TAPER  TIME  SERIES 


PLOT  TIME  SERIES 


PERFORM  FFT  OF  TIME  SERIES 


PLOT  SIGNAL 
SPECTRUM 


CORRECT  SIGNAL  SPECTRUM  FOR 
INSTRUMENT  RESPONSE 


PLOT  CORRECTED 
SIGNAL  SPECTRUM 


CALCULATE  NOISE  AND  SIGNAL 
POWER  SPECTRA:  SIGNAL-TO- 
NOISE  RATIO 

LOOP  OVER  DIFFERENT 


LTER  FREQUENCIES 


PLOT  NOISE  AND 
SIGNAL  POWER; 
SIGNAL-TO-NCISE 
RATIO 


NARROW  BAND  FILTER  (NBF) 
CORRECTED  SIGNAL  SPECTRUM 


CALCULATE  QUADRATURE 
SPECTRUM  BY  HILBERT  TRANSFORM 


Figure  la.  MARS  flowchart. 
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Figure  la.  (continued) 
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LOOP  OVER 
SEQUENCE  OF 
SEISMOGRAMS 


READ  AND  EDIT  INPUT  DATA 


Figure  lb.  Flowchart  indicating  principal  mathematical 
operations  in  the  MARS  program. 
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Figure  lb.  (continued) 
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Figure  lb.  (continued) 
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COMPUTE  POLARIZATION  FILTER  FACTORS  FOR  EACH  PAIR 

OF  COMPONENTS; 

E.G.,  THE  VERTICAL-RADIAL  (Z-R) 
POLARIZATION  FACTORS  ARE 


p£k)  (t)  -  sin  A  $£k)  (t) 


CIRCULAR  FUNCTION 


(t)  -  cos  A  $j;k)  (t) 


LINEAR  FUNCTION 


SR(k) (t)  -  A^k) (t)/A^k5  (t) 


ELLIPSE  RATIO 


ALF°°(t)  =  arctan  [A^k)  (t) /A^k)  (t)  ]  APPARENT  P- WAVE 


EMERGENCE  ANGLE 


WHERE  A$^k>  »  $^k)  (t)  -  (t) 


COMPUTE  POLARIZATION  AND  DIRECTIONAL  FILTER  DETECTORS 
FOR  SPECIFIC  TYPES  OF  WAVE  MOTION; 

E.G.,  THE  P-WAVE  DETECTOR  IS 

P  00  (t)  *  cos**1  A  5^k)  (t) 

•  cos  2  [arctan  (A^k'/A^‘<J)  -  a,  ]  , 
if  Pjfk)  (t)  >  0 

WHERE  M.  AND  M,  ARE  INTEGERS  AND  a.  IS  THE  APPARENT 

EMERGENCE  ANGLE. 


(continued) 
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Figure  lb.  (continued) 
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Seismic  data  are  read  into  the  program  in  the  form  of 
time  series  (Figures  la  and  lb)  generally  of  about  500  to  2000 
points  in  length.  The  data  are  then  optionally  detrended, 
demeaned  and  tapered.  The  program  then  selects  the  smallest 
power  of  two  which  is  greater  than  the  number  of  points  input 
and  performs  a  discrete  Fourier  transform  using  the  Cooley- 
Tukey  algorithm.  The  signal  spectrum  is  corrected  for  the 
appropriate  instrument  response;  i.e.,  the  signal  transform 
is  divided  by  the  instrument  transfer  function. 


Referring  to  Figures  la  and  lb,  the  signal  is  next 
filtered  in  the  frequency  domain  by  multiplication  by  a  narrow 
band  Gaussian-shaped  filter.  This  particular  filter  form  is 
selected  to  satisfy  two  goals:  (1)  minimum  width  in  the  fre¬ 
quency  domain,  and  (2)  maximum  ripple  suppression  in  the  time 
domain.  As  a  consequence  of  the  sampling  theorem,  or  uncer¬ 
tainty  principle,  one  cannot  simultaneously  satisfy  these  two 
goals  to  arbitrary  precision.  The  filter  employed  was  selected 
for  its  optimal  time  and  frequency  domain  characteristics 
within  this  basic  limitation.  Once  the  signal  has  been  narrow 
band  filtered,  the  resulting  complex  spectrum  is  then  inverse 
Fourier  transformed  into  the  time  domain,  to  produce  what  will 
hereafter  be  referred  to  as  the  filtered  signal. 

The  narrow  band  filtered  signal  will  appear  as  a  quasi- 
sinusoiaal  carrier  wave  contained  within  a  smooth  envelope. 

The  next  step  in  the  program  (Figures  la  and  lb)  is  to  con¬ 
struct  the  envelope  function  by  means  of  the  Hilbert  trans¬ 
form.  In  particular,  a  quadrature  signal  is  formed  by 
multiplying  the  transform  of  the  filtered  signal  by  -i  sgn  (oj) 
and  then  inverting  this  to  the  time  domain  by  an  inverse  trans¬ 
formation.  The  envelope  function  is  then  constructed  by  taking 
the  square  root  of  the  sura  of  the  squares  of  the  filtered  sig¬ 
nal  and  its  quadrature.  The  maximum  of  the  envelope  function 
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occurs  at  the  time  of  arrival  of  energy  at  the  center  frequency 
of  the  filter  and  the  amplitude  of  the  maximum  is  proportional 
to  the  spectral  amplitude  of  the  filtered  signal  at  the  center 
frequency  of  the  filter.  The  narrow  band  filtering  procedure 
can  be  performed  on  a  particular  component  seismogram  at  a 
number  of  different  frequencies  within  some  band  of  interest. 
Correlation  of  the  resulting  envelope  functions  indicates  the 
arrival  times  of  the  various  frequency  components. 

The  instantaneous  frequency  and  phase  are  computed  from 
the  filtered  signal  and  its  quadrature  signal;  both  signals 
are  stored  for  subsequent  use  in  polarization  filtering  with 
additional  components  of  ground  motion  when  available. 

The  output  of  the  narrow  band  filter  processing  is  a 
table  including  the  envelope  peak  amplitudes  (A^)  and  asso¬ 
ciated  arrival  times  (t^)  for  each  center  frequency  (f^) . 
Generally  there  are  many  A^,  T  pairs  for  each  filter, 
depending  on  the  complexity  of  the  signal. 


(2) 


a  *  1  _k 

f  ~  >/En2  Q 


which  means  that  68.3  percent  of  the  area  under  the  filter  is 
contained  in  the  frequency  band  f^  -  ac  <_  f^  <  f^  +  crf. 

It  is  easily  shown  that  the  time  domain  expression  of 
the  filter  (Equation  1)  is 


*2t2 


F(k)  (t)  =  —  cos  (27rfkt) 


(3) 


The  Gaussian  envelope  function  then  has  standard  deviation 


_  y/£n2  _Q_ 
T  “  2 it  f, 

Jv 


(4) 


We  see  that  the  width  of  the  narrow  band  filter  is  con¬ 
trolled  by  the  parameter  Q  which  appears  in  the  definition  of 
a.  From  Equations  (2)  and  (4)  it  is  clear  that  the  narrower 
the  filter  in  frequency,  the  wider  in  time  and  vice  versa. 

For  high  Q  we  get  an  excellent  estimate  of  the  spectral  ampli¬ 
tude  for  a  broad  time  window;  there  will  be  few  A^,  t^  pairs 
for  each  filter.  For  low  Q  the  spectral  estimate  is  less  pre¬ 
cise  but  the  tg  are  more  accurate  and  there  are  generally  many 
peaks  for  each  filter. 

A  choice  of  a  preferred  Q  was  made  by  trial  and  error 
testing  of  different  values  on  both  synthetic  and  observed 
seismograms.  For  short  period  data  sampled  20  times/second, 
we  have  adopted  the  following  Q  values  for  all  the  seismograms 
processed  in  this  experiment: 
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fk  1  0.70 


•Q  -  12.5  fk, 

Q  -  8.00,  fk  <  0.70.  (5) 

For  fk  >_  0.70  we  compute  from  Equation  (2)  that 
af  =  0.096  Hz. 

Thus,  the  width  of  the  filter  within  one  standard  deviation  is 
0.192  Hz.  For  fk  <  0.70  Equation  (2)  gives  cf  *  0.15  fk  and 
the  filter  is  narrower. 

From  Equation  (4)  we  have  for  fk  ^  0.70  that 
aT  *  1.66  seconds. 

The  filter  width  is  therefore  approximately  three  seconds  for 
these  fk,  using  standard  deviation  as  the  criterion,  and  some¬ 
what  wider  at  lower  center  frequencies. 

Figure  2  shows  the  outputs  of  a  suite  of  narrow  band 
filters  applied  to  two  event  seismograms ,  event  47  in  Figure 
2a  and  event  21  in  Figure  2b.  As  shown  later,  event  47  is 
identified  as  a  Kurile  earthquake  while  event  21  is  identified 
as  a  presumed  explosion  north  of  Lake  Baikal.  The  original 
unfiltered  seismograms  are  the  topmost  traces  in  each  figure. 
The  center  frequencies  of  the  narrow  band  filters  range  from 
4.0  to  0.5  Hz.  In  terms  of  sigr.al-to-noise  ratio  (S/N),  the 
basic  features  of  the  VFM  discriminant  can  be  seen  here.  For 
instance,  the  presumed  explosion  signal  in  Figure  2b  is  quite 
evident  over  the  frequency  range  3.5  to  0.8  Hz.  The  earthquake 
signal  (Figure  2a)  does  not  persist  above  the  ambient  noise 
level  at  frequencies  higher  than  2.75  Hz,  whereas  it  exhibits 
a  higher  S/N  at  1.0  Hz  and  0.8  Hz  than  the  presumed  explosion. 

A  final  point  to  be  noted  frjm  Figures  2a  and  2b  i3  the 
approximate  time  alignment  of  the  maximum  narrow  band  filter 
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Figure  2a.  Examples  of  the  outputs  of  a  suite  of  4.0  to  0.5 
Hz  narrow  band  filters  applied  to  a  seismogram 
(shown  at  top  of  figure)  from  a  Kuril  earthquake 
recorded  at  a  station  at  teleseismic  distances. 
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Figure  2b.  Examples  of  the  outputs  of  a  suite  of  narrow 
band  filters  applied  to  a  seismogram  (top  of 
figure)  of  a  presumed  explosion  recorded  at 
the  same  station  as  in  Figure  2a. 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


amplitudes,  with  high  S/N,  corresponding  to  the  signal  arrival. 
This  observation  forms  the  basis  for  the  automated  detection- 
discrimination  algorithm  described  in  the  following  subsection 
of  this  report. 


2.1.2  The  Automated  Detection  Algorithm 

The  decomposition  of  a  signal  wave  train  into  quasi¬ 
harmonic  signals  by  the  application  of  phaseless  Gaussian- 
shaped  narrow  band  filters  provides  the  means  of  determining 
arrival  time,  amplitude  and  phase,  all  as  functions  of  fre¬ 
quency.  This  is  the  basic  signal  information  that  can  be 
used  to  detect  a  given  type  of  signal  in  terms  of  its  disper¬ 
sion  characteristics  and  to  obtain  its  spectrum  as  well  as 
its  time  and  amplitude  relationship  with  respect  to  other 
signals  present  in  a  complex  wave  train. 

The  basic  approach  used  in  the  MARS  program  is  to 

identify  patterns  from  the  signal  and  noise  information  as  it 

is  expressed  in  the  t  -f  olane  (the  group  arrival  time,  t  , 

g  g 

versus  frequency  plane) .  The  pattern  to  be  searched  for  in 

the  t^-f  plane  corresponds,  in  the  case  of  a  body  wave,  to  a 
(nearly)  undispersed  signal,  with  spectral  amplitude  signifi¬ 
cantly  above  background  in  a  frequency  range  corresponding  to 
some  fraction  of  the  total  band.  This  frequency  band  will  be 
in  a  range  where  the  signal  power  is  expected  to  be  highest 
relative  to  noise.  (This  means  we  make  use  of  the  matched 
filtering  concept,  in  the  sense  that  we  know  roughly  what 
spectral  content  we  expect  for  the  signal.  This  concept  is 
also  used  when  we  look  for  undispersed  signals  or  signals  of 
known  dispersion  characteristics.)  Thus,  in  the  t  -f  plane 
we  search  for  a  straight,  horizontal  line-up  of  envelope 
maxima  in  a  selected  frequency  band,  using  the  largest  maxima 
as  the  beginning  point  in  such  a  pattern  recognition  procedure. 
This  implicitly  uses  a  threshold  detection  criteria  in  a 
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detection  band,  since  by  starting  with  the  largest  envelope 
maxima  we  are  essentially  requiring  that  the  signal  power  be 
above  the  background  level  in  at  least  a  part  of  the  detec¬ 
tion  band. 

Basically,  then,  by  looking  in  this  t  -f-A  -$  space, 
we  apply  criteria  based  on  properties  of  the  expected  signal, 
namely  its  expected  dispersion  and  spectral  content,  in  order 
to  recognize  a  signal  pattern  and  thereby  to  detect  the  sig¬ 
nal.  An  example  of  how  this  is  accomplished  is  shown  in 
Figure  3.  In  this  figure  we  plot  the  times  (tg)  of  the 
envelope  maxima  from  narrow  band  filter  outputs,  with  of  the 
order  of  N  =  20  filters  used  so  that  the  signal  frequency 
content  is  sampled  at  about  20  points,  f  .  With  each  envelope 
maxima  point  in  the  plane,  there  is  also  an  associated  (spec¬ 
tral)  anrolitude  A  (f  ) ,  instantaneous  phase  $  (f  )  and  an 
*  g  n  c  g  n 

instantaneous  frequency  (dsjj/dt)^  .  Thus  the  t  -f  plot  corre- 

wg  g 

sponds  to  a  multidimensional  display  of  spectral  content  and 
energy  arrival  time  for  a  given  segment  of  a  time  series. 
Normally  either  1024  or  2048  points  are  used  for  time  segments, 
and  for  the  short  period  seismic  data  in  the  discrimination 
data  base  this  corresponded  to  a  50  to  100  second  segment 
which  was  processed  in  each  pass.  For  on-line  continuous 
processing,  overlapping  time  segments  would  be  used. 

A  sub-band  within  the  entire  frequency  range  covered 
by  the  set  of  filters  is  shown  in  Figure  3  and  is  used  as  a 
"detection  band,"  that  is  a  frequency  band  within  which  a 
signal  pattern  (straight  horizontal  line  locus  of  envelope 
maxima  in  the  case  of  an  undispersed  body  wave)  is  sought. 

This  band,  from  f*  to  f*,  is  selected  externally,  based  on 
the  expected  signal  frequency  character.  The  largest  envelope 
maxima  within  this  band  are  flagged  (denoted  by  0  in  the 
figure)  and  used  to  compute  a  mean  "signal"  arrival  time 
(i.e.,  group  time)  t”  for  the  maximum  power  arriving  in  this 
frequency  range.  An  acceptance  window  in  time,  for  which 
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Figure  3. 


Typical  tg-f  plane  representation  of  a  time  series  seg¬ 
ment  (0  to  36  sec)  when  the  signal- to-noise  ratio  is  low. 
Only  the  largest  envelope  maxima  (X)  and  second  largest 
maxima  (O)  are  3hown.  Other,  numerous,  smaller  envelope 
maxima  are  normally  scattered  throughout  the  t„-f  plane, 
but  for  this  illustration  they  have  been  omitted.  Each 
filter  "output  line"  normally  displays  about  ten  such 
peaks. 
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maxima-  can  actually  be  associated  with  an  undispersed  signal, 

+  —  + 

is  constructed  using  the  relation:  t*  ■  t+  t"  +  aAt;  where  At 

w  g  —  cr  — 

is  the  time  uncertainty  associated  with  the  envelope  maximum 
time  t  for  a  filter  output  at  center  frequency  f  and  half 
power  band  width  Af.  In  particular,  AtAf  >_  1/4II  is  the  theo¬ 
retical  uncertainty  relation  and  since  Af/f  =  Q-^,  then 
At  >_  1/411  (Q/f).  Further  we  use  t*  to  denote  the  standard 
deviation  in  the  time  data  used  to  compute  t  and  make  use  of 
it  to  define  the  acceptance  time  window  about  t  as  well. 

Hence,  taking  an  appropriately  chosen  constant  a  near  unity, 
then  two  time  window  boundary  lines  can  be  defined  to  give 
the  time  window;  that  is:  t_  +  t_  +  (a/43I)  (Q/f)  and 
t  -  t^  -  (a/4H) (Q/f).  All 'of  the  largest  maxima  within  such 
a  window  are  then  taken  as  acceptable  undispersed  "signal" 
group  arrivals  and  those  outside  are  rejected. 

The  largest  maxima  within  the  initially  defined  window 
are  now  used  to  recompute  a  new  group  arrival  time,  t  ,  and 
signal  acceptance  window.  Finally,  using  all  those  filter 
peaks  that  lie  within  the  tg  acceptance  window  and  exceed  the 
ambient  noise  level,  we  define  a  final  group  arrival  time,  t  . 

g 

The  t_  estimate  corresoonds  to  the  arrival  time  of 

g 

maximum  signal  energy.  For  purposes  of  multiple  phase  detec¬ 
tion  and  event  locations  we  added  an  algorithm  to  MARS  that 
gives  estimates  of  phase  arrival  times.  This  is  accomplished 
by  flagging  the  earliest  filter  peaks  (envelope  maxima)  that 
occur  within  the  final  signal  acceptance  window.  There  must 
be  a  minimum  number  of  such  peaks,  typically  three,  all  of 
which  exceed  the  mean  noise  level  at  the  corresponding  fre¬ 
quencies  by  at  least  one  standard  deviation. 

Figures  4  through  7  are  a  sequence  of  examples  of  the 
application  of  the  automated  detection  and  first  arrival 
algorithms  to  seismograms  for  two  of  the  events  in  the  dis¬ 
crimination  data  base;  events  47  and  81.  The  events  are 
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Figure  5.  Seismogram  (bottom)  written  at  TataJina,  Alaska,  for  the  same 
presumed  explosion  as  in  Figure  4. 
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igure  6.  Bottom  portion  is  a  seismogram  recorded  at  a  station 
in  Bluff,  Alaska,  from  an  event  in  the  Kurils.  The 
top  half  is  a  plot  of  filter  center  frequency  versus 
group  arrival  time  (sec).  The  t|  with  arrow  attached 
shows  where  the  automated  detection  and  first  arrival 
algorithms  picked  the  signal  arrival  time. 
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Bottom  portion  is  a  seismogram  recorded  at  NORSAR 
for  the  same  Kuril  earthquake  shown  in  Figure  6. 
The  top  portion  of  this  figure  is  the  fc  versus  t 
plot  determined  for  this  time  series  by  MARS.  " 
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shown  in  the  order  of  decreasing  signal-to-noise  ratio.  The 
bottom  portions  of  these  figures  show  the  seismograms  that  were 
processed  with  MARS.  In  each  case,  the  recordings  are  short- 
period  vertical  component  seismograms.  The  corresponding  digital 
time  series  were  100  seconds  in  duration  with  at  least  50  seconds 
of  background  noise  preceding  the  signal  arrival. 

The  top  portions  of  Figures  4  through  7  are  the  final 

t  -f  planes  that  are  arrived  at  after  iterating  on  the  filter 
g  c 

peaks  as  described  above.  The  group  arrival  time,  measured 
from  the  event  origin  time,  is  plotted  along  the  abscissa  and 
filter  center  frequency  along  the  ordinate  in  each  of  these 
figures.  The  relative  amplitudes  of  the  filter  envelope  maxima 
are  represented  by  the  numbers  0  to  9  where,  for  instance,  a  7 
corresponds  to  an  envelope  peak  that  is  between  70  and  80  per¬ 
cent  of  the  largest  envelope  peak,  designated  by  an  (*) ,  at  the 
same  frequency.  The  +  and  -'s  define  the  signal  acceptance 
window  and  give  a  measure  of  the  uncertainty  in  time  resolution 
due  to  scatter  {in  time)  of  the  envelope  peaks  and  the  inher¬ 
ent  limitation  associated  with  the  narrow  band  filtering 
procedure. 


For  most  of  the  event  seismograms  (Figures  4,  5  and 
6)  the  signal-to-noise  ratio  (S/N)  is  sufficiently  high  so 
that  it  is  obvious  from  an  examination  of  the  figures  that  the 
arrival  times  picked  by  the  algorithm  correspond  to  those  that 
an  analyst  would  pick.  Figure  7,  however,  is  an  entirely  dif¬ 
ferent  matter  and  the  close  agreement  between  t|  and  the  pre¬ 
dicted  arrival  time  (tp^,  based  on  the  known  event  location, 
origin  time  and  Herrin  (1968)  travel-time  tables)  is  indica¬ 
tive  of  the  power  of  the  narrow  band  filter  approach  to  signal 
detection. 


An  example  of  the  multiple  phase  recognition  capabilities 
of  the  automated  detection  algorithm  implemented  in  MARS  is 
shown  in  Figure  8.  The  bottom  portion  of  this  figure  is  a 
short-period  vertical  component  seismogram  written  at  the 
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Figure  8.  Example  of  timing  of  multiple  phase  arrivals  by  MARS.  The  seismogram 
the  bottom  is  for  the  NTS  event,  Camembert,  written  at  Houlton,  Maine 
The  top  is  the  t  -f  plot  for  a  series  of  narrow  band  filters. 


station  in  Houlton,  Maine  (HNME)  for  the  NTS  explosion 

Camembert.  The  three  arrows  point  to  the  arrival  times  of 

three  separate  phases  picked  by  MARS.  The  phase  arrival 

times  are  based  on  the  t  -f  plot  shown  in  the  top-half  of 

<3 

this  figure. 


Table  1  compares  the  lag  times  picked  by  MARS  for 
Camembert  and  Mast  with  those  used  in  synthetic  seismogram 
calculations  for  these  particular  events.  The  synthetic 
seismograms,  which  matched  the  observed  very  closely  (Bache, 
1978)  required  the  superposition  of  three  separate  phases 
corresponding  to  the  direct  P-wave,  a  second  phase  with  a 
delay  appropriate  for  pP  and  a  third  phase  that  could  be 
either  associated  with  spallation  or  tectonic  stress  release. 
Whatever  the  correct  explanation  for  the  later  arrivals  here, 
the  remarkable  agreement  between  the  observed  and  calculated 
lag  times  implies  that  we  are  well  on  our  way  toward  developing 
an  automatic  multiple  phase  identification  technique. 


Once  a  final  signal  acceptance  window  has  been  defined 
(e.g.,  the  +  and  -'s  in  the  tg-f  planes  in  Figures  4  through 
8),  the  next  step  is  to  obtain  estimates  of  the  signal  ampli¬ 
tudes  at  the  different  frequencies  for  subsequent  magnitude 
calculations.  To  do  this  a  combined  amplitude-time  scoring 
scheme  is  applied  to  each  of  the  filter  peaks  that  occur 
within  the  signal  acceptance  window.  The  time  scoring  is  done 
with  respect  to  the  t  estimate  indicated,  for  example,  by  the 
B's  at  the  top  and  bottom  of  each  of  the  t  -f  signal  acceptance 
windows  in  Figures  4  through  8.  The  amplitude  scoring  is 
straightforward  with  the  largest  filter  peak  receiving  a  nor¬ 
malized  amplitude  score  of  one.  At  a  particular  frequency  the 
highest  total  score  will  be  awarded  the  largest  peak  that 
occurs  at  time  t^.  With  the  signal  amplitudes  defined  over 
the  frequency  band  of  interest  in  this  experiment  (i.e.,  0.4  Ez 
to  5.0  Hz  for  the  20  samples/sec  data  and  0.3  to  3.5  Hz  for  the 
10  samples/sec  data) ,  we  turn  to  the  discrimination  algorithm. 
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TABLE  1 

COMPARISON  OF  PHASE  LAG  TIMES 


MAST 

Synthetic 
Observed  (MARS) 


P2  -  P, 
(sec) 


P3  -  P1 
(sec) 


0. 67  1.35 

0.69  1.25 


CAMEMBERT 

Synthetic  0.64  1.92 

Observed  (MARS)  0.64  1.93 
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2.1.3  The  Discrimination  Algorithm 

The  results  of  one  of  our  earliest  (Savino  and  Archambeau, 
1974)  applications  of  the  variable  frequency  magnitude,  or  ra^(f) , 
discriminant  to  seismic  events  are  shown  in  Figure  9.  The  data 
shown  are  only  two  of  the  frequency  dependent  magnitudes  that 
characterize  the  events,  but  for  the  teleseismic  distances 
(A  =  60  to  90°)  and  the  type  of  receiver  (LASA  -  full  array  beam) 
the  two  frequency  dependent  magnitudes  at  0.45  and  2.25  Hz  pro¬ 
vide  the  best  separation  of  the  explosion-earthquake  populations 
in  the  C f )  parameter  space. 

An  important  point  about  the  results  in  Figure  9  con¬ 
cerns  the  influence  of  background  seismic  noise  on  the  event 
populations.  Some  of  the  scatter  in  both  the  earthquake  and 
explosion  populations  is  due  to  noise.  This,  of  course,  is 
relatively  more  important  for  events  of  small  magnitude,  since 
the  signal  power  is  low  relative  to  the  noise  power.  An  addi¬ 
tional  effect  is  mixing  of  the  populations  at  the  very  low 
magnitudes. 

In  view  of  the  desirability  of  estimating  noise  prop¬ 
erties  and  then  correcting  for  noise  contamination,  the  signal 
analysis  program,  as  originally  written,  obtained  noise  spectral 
estimates  and  used  these  to  obtain  estimates  of  the  signal  spec¬ 
tral  magnitudes.  (In  the  original  case,  the  average  noise  spec¬ 
tral  amplitude  at  the  appropriate  frequency  was  subtracted  from 
the  signal  spectral  amplitude.)  The  nature  of  the  noise  has 
been  investigated  for  each  of  the  events  making  up  the  event 
set  shown  in  Figure  9  using  the  signal  analysis  program,  and 
it  was  found  that  noise  time  series,  when  viewed  in  an  m^Cf) 
plane  such  as  that  in  Figure  9,  occupies  the  same  region  as 
do  the  small  magnitude  earthquakes .  This  means  that  when  the 
noise  contaminates  a  small  explosion  signal,  the  result  will 
be  such  as  to  make  the  explosion  appear  to  be  more  earthquake¬ 
like  and,  for  low  signal-to-noise  ratios,  to  actually  move  the 
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Figure  9.  Spectral  magnitudes,  mb(f),  computed  at  0.45  Hz  and 
2.25  Hz.  The  presumed  explosions  numbered  1  and  2 
occurred  at  Novaya  Zemlya,  all  the  other  explosions 
occurred  in  eastern  Kazakhstan. 
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explos-ion  point  into,  or  very  near,  the  earthquake  population. 
Further,  the  noise  contamination  will  move  the  earthquake 
points,  but  always  within  the  earthquake  population  and  not 
outside  of  it. 

This  is  illustrated  in  Figure  10a  where  a  number  of  the 
small  events  observed  at  LASA  and  shown  in  Figure  9  have  been 
detected  at  a  Norway  array  and  processed  and  plotted  in  the 
m^ff)  plane  for  f  =  0.6  Hz  and  f  =  5.0  Hz.  Again,  in  this 
plane,  we  observe  convergence  of  the  populations.  The  arrows 
on  the  smallest  explosion  data  points  indicate  the  direction 
in  which  a  noise  correction  will  move  these  points  in  this 
magnitude  plane.  Figure  10b  shows  the  effect  of  the  noise 
correction  applied  to  all  the  events,  including  earthquakes 
as  well  as  explosions.  (In  this  figure  only  the  log  of  the 
amplitude  is  plotted,  the  m^  plot  would  be  identical  except 
for  a  scale  change.)  It  is  evident  that  the  explosion  popula¬ 
tion  is  now  well  separated  from  the  earthquake  population.  The 
arrows  on  the  two  small  explosions,  both  with  regular  m^  mag¬ 
nitudes  somewhat  less  than  4.0,  indicate  that  it  is  likely  that 
the  estimate  of  the  noise  contamination  is  lower  than  that 
actually  present  and  that  the  true  event  points  are  still  lower 
in  the  direction  indicated. 

While- the  application  of  a  noise  "correction"  to  the 
data  in  Figure  10  was  observed  to  improve  discrimination,  the 
form  of  this  correction,  namely  the  subtraction  of  an  average 
noise  level,  was  known  to  produce  a  biased  (high)  estimate  of 
the  probable  noise  level  occurring  during  the  signal  window 
and  also  did  not  properly  account  for  the  phase  of  the  noise 
relative  to  that  cf  the  signal.  Thus,  a  more  consistent  and 
logical  algorithm  for  estimating  the  influence  of  seismic 
noise  on  a  transient  signal  was  formulated  and  implemented  in 
MARS. 

In  principle,  the  noise  correction,  as  now  applied,  con¬ 
sists  of  a  deterministic  component  and  a  statistical  component. 
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filter  amplitudes  with  noise  correction 
same  event  population  plotted  in  Figure 
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The  deterministic  component  is  based  on  a  superposed  pulse 
model  for  the  noise#  where  noise  is  defined  as  all  energy  or 
group  arrivals  not  identified  with  the  particular  signal  being 
considered.  Given  this  definition,  the  "noise"  can  be  made  up 
of  what  would  ordinarily  be  considered  signal  (e.g.,  the  coda 
of  the  first  arrival  P-wave) ,  as  well  as  normal  background 
seismic  noise.  The  particular  form  of  the  correction  that 
treats  this  "local  seismic  noise"  is  termed  deterministic 
because  the  effects  of  both  the  amplitude  and  instantaneous 
phase  of  the  "noise"  on  the  signal  can  be  calculated,  at  least 
to  first  order. 

As  described  earlier,  the  output  of  the  narrow  band 

filtering  process  consists  of  maxima  of  envelope  functions 

as  a  function  time.  The  deterministic  noise  correction  is 

formulated  as  follows.  Let  A*(f)  be  the  measured  envelope 

amolitude  associated  with  a  signal  of  interest,  and  t*(f)  be 

5 

the  energy,  or  group,  arrival  time.  In  addition,  let  (A  (f )  } 

N  n 

be  a  set  of  noise  peaks  with  group  times  (tg(f)}  such  that, 
either:  t*  -  ot  <  t**(f)  <  t*  -  it,  or  t*  +  it  <  tN(f)  <  t*  +  <$t; 

g  _g_g  g  _g  _g 

where: 


=  10  &n2.  „  =  _1_ 

and  suppose  that  there  are  M  such  noise  peaks.  Here  Au  is  the 
half  power  band  width  of  the  Gaussian  filters  that  are  used  in 
MARS. 

We  compute. the  "local  noise  corrected  signal  amplitude," 
A**  ( f ) ,  from: 


37 


SYSTEMS.  SCIENCE  AND  SOFTWARE 


A*  cos 
9  9  9 


A** (f )  = 


M 

E 

m=l 


-  2  exp 


~Ah>,  /+■*  -  tm)  ^ 

4S  ltg  V 


cos  (t*)  ) 

P  9  j 


M 


A*  sin  4<*(t*)  -  52  ^ 

g  p  g  nfci  ^ 


(m) 


exp 


.  (t* 
43  ltg 


tgm)  ) 


sin  f[m) (t*)| 

p  g  1 


1/2 


where: 


3  = 


In  2 


TJ(t5) 


A*(t*> 
9  9 


^(ra)  (ta} 

p  g 


,?(m)  (t*} 

p  g 


AJJm)  (f) 


Instantaneous  phase  of  the  signal  at 
the  envelope  peak  time  (group  time) 

t|(f)- 

Amplitude  of  the  envelope  at  the  peak 
in  the  envelope  function,  occurring  at 

t|(f). 

f  K 

Instantaneous  phase  of  the  m  *  "noise" 
pulse  at  the  envelope  peak  time  (group 
time),  t®(f). 

u>(t*  -  t“)  +  (t®)  =  the  noise 

o  g  g  P  9 
phase  at  t|,  the  signal  group  time. 

Amplitude  at  the  envelope  peak  for  the 
mfc^  noise  pulse  (occurring  at  the  time 

t”(f>>. 
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Here: 


M 


sn  *  2  j4ml  «p 


z&a i .  (t*  -  tm)2 

46  v 


cos  (t*) 
p  g 


is  the  deterministic  noise  correction  for  the  signal  with  group 
time  t*,  while 


2 


sin  (t*) 

p  g 


is  the  quadrature  component  corresponding  to  oA^.  We  then 
compute  both  of  these  and  evaluate  the  (deterministic)  noise 
correction  J  [  as: 


2  .  2  1/2 
ISAj,!  »  [{«J V  +  (SAj,}  1 

and  save  for  later  use  in  describing  the  noise  population  in 
the  ra^f)  plane.  This  correction,  while  deterministic,  is 
only  to  first  order. 

Associated  with  A|*(f)  will  be  an  uncertainty  due  to 
the  impossibility  of  resolving  and  correcting  for  noise  peaks 
within  the  range  At  »  l/2Au  on  either  side  of  the  signal  group 
time  (t*) ,  at  a  frequency  f.  In  addition,  we  will  want  to 
include  in  the  uncertainty  attached  to  the  estimate  A**,  the 

g 
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impossibility  of  making  an  exact  deterministic  correction  for 
noise  contamination  from  the  "local"  noise  pulses.  This  uncer¬ 
tainty  will  be  +  AAjj,  with  AA^j  given  by: 

-  <S>  l£>  V£>  =  A  <£>  *8 

o  o 

where 

Tq  =  Total  (standard)  time  window  being  processed 
(typically  above  100  sec) . 

L  «  Number  of  noise  peaks  in  Tq.  (Use  the  pre¬ 
viously  analyzed  time  window  where  all 
noise  peaks  have  been  identified. ) 

A^(f)  =  Mean  of  the  envelope  amplitudes  at  the  noise 
peaks  in  the  window  Tq  ,  i . e . , 

L 

V*)  *  2  AjS1’  (f)/t  . 

N  1=1 

(Note  that  L  a  iu  Tfl  as  iu  •*  0,  so  that 
AAjj  •+■  Ajj  as  Aw  -*■  0.  Further,  *  0  as 
Aw  •*  0  by  definition  of  the  time  interval 
for  the  oA^  correction.  Therefore  as  iu  +  0, 
and  the  filter  Q  become  infinite,  we  get  the 
usual  noise  correction  of  a  Fourier  spectrum.) 

Thus  the  signal  spectral  amplitude  will  be 
described  by  A|  (f)  +  AA^(f)  . 

The  final  step  in  this  formulation  is  the  computation 
of  the  noise  corrected  instantaneous  phase.  This  is  given  bv: 
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y**(t*(f) ) 
P  9 


+  tan 


-1 


y*  ~ 

y*  -  <5^ 


and  SA^  are  as  defined  above  while 
y*  =  A*  cos  V*(t*) 

g  p  g 

y*  =  +  A*  sin  ¥*(1*) 

where  y*  is  the  quadrature  signal. 

The  application  of  the  noise  correction  derived  above 
is  based  on  a  magnitude  relationship  similar  to  the  one 
originally  proposed  by  Gutenberg  and  Richter  (1956): 

m^f)  »  log^Q  [Af 3  +  b 

where  b  is  the  distance  correction  factor.  Using  this  rela¬ 
tionship  we  compute  m^Cf)  values  for  the  signal  of  interest 
from  A*  ( f )  ,  A**(f),  A*Mf)  +  AA^  and  A**(f)  -  AA^.  The  m^f) 
are  computed  at  frequencies  corresponding  to  the  center  fre¬ 
quencies  of  the  entire  set  of  narrow  band  filters  being  used 
(i.e.,  typically  40  filters  covering  the  band  0.4  to  5  Hz). 

For  discrimination  purposes,  however,  two  sub-bands  are 
defined:  a  low  frequency  set  { f^}  and  a  high  frequency  set 

{ f H  > -  The  high  and  low  frequency  bands  are  defined  in  Table 
2.  In  order  to  obtain  more  stable  magnitude  estimates 
than  those  based  on  individual  filter  center  frequencies, 
fitted  values  at  a  single  fixed  high  (fH)  and  a  fixed  low 
(f^)  frequency  are  used.  For  the  high  frequency  band,  the 
least  squares  polynomial  of  degree  two  which  best  approximates 

the  weighted  set  of  log. „  A**(f)  values  over  the  band  { f „ >  is 

1U  y  tt 

first  computed.  The  weight  factors  w,  (f)  which  measure  the 


where  <5A^ 
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TABLE  2 

HIGH  AND  LOW  FREQUENCY  BANDS 


<v 

ffL> 

Sampling 

Rate 

Number  of 
Filters 

5.0  -  2.0 

1.4  -  0.40 

20 

40+ 

3.5  -  2.0 

1.4  -  0.30 

10 

40  tt 

f  Twenty 

samples/second 

stations 

tt  NORSAR  and  LASA 
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signal* information  content  or  quality  of  the  signal  information 
at  the  frequency  f^  are  given  by 

wk(f)  =  [A*  { f )  /Ajj  ( f )  ]  2  . 

The  fitted  value  at  the  single  high  fixed  frequency  fg  is  then 
determined  and  used  to  derive  a  magnitude  estimate.  A  similar 
procedure  is  applied  to  derive  the  low  frequency  estimate. 

The  magnitude  relationship  is  given  by 


mb(F)  =  log1Q  [A(F>  •  F]  +  b  , 


where  F  is  the  single  fixed  frequency  f„  or  fL  and  A (F)  is  the 
fitted  value  computed  from  the  weighted  least-squares  polynomial 
defined  above  for  the  set  of  noise  corrected  amplitudes  A**(f). 

Finally,  the  uncertainty  in  the  magnitudes  at  f^  and  fg 
due  to  the  uncertainty  in  the  noise  correction  is  also  computed. 
These  are  given  by: 

Am£(F)  =  log1Q  C(A(F)  +  AA^)  •  F] 

+  b  -  m^ (F) , 

An£(F)  *  log1Q  C(A(F)  -  AA^)  •  F] 

+  b  -  nih  (F)  , 

where  A (F)  is  the  fitted  value  for  the  noise  corrected  ampli¬ 
tude  A**  at  frequency  F,  AA^  is  the  noise  uncertainty  at  F, 
and  F  =»  fH  or  f^.  Note  that  Am£  ^  - 

none  of  the  Am,  values  obtained  for  the  two  sets  of  discrimina- 
o 

tion  frequencies  will  be  equal  in  general.  We  get  from  this 
then  four  distinct  Am^  values. 
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The  results  of  applying  the  newly  formulated  noise  cor¬ 
rection  to  event  m^ff)  data  are  summarized  in  Figures  11a  and 
lib.  Figure  11a  shows  the  behavior  of  typical  event  (uncor¬ 
rected  for  noise)  and  noise  populations  in  the  m^Cf)  plane. 
This  figure  is  a  generalization  of  earlier  results  previously 
reported  (Savino,  et  al. ,  1975;  Rodi,  et  al. .  1978).  Figure 
lib  demonstrates  the  manner  in  which  the  earthquake  and  ex¬ 
plosion  populations  are  expected  to  separate  when  the  deter¬ 
ministic  and  statistical  noise  corrections,  together  with  the 
uncertainty  inherent  in  both  these  corrections,  are  applied. 
The  enhanced  separation  of  populations  is  especially  signifi¬ 
cant  in  the  low  m^(f)  range  where  noise  plays  an  important 
role.  This  is  the  procedure  that  was  routinely  used  in  the 
discrimination  experiment  for  the  computation  of  the  variable 
frequency  magnitudes. 
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♦  ■  ■ 


NOISE  POPULATION  BOUNDARY 


HIGH  FREQUENCY  BODY  WAVE  MAGNITUDE  [m^ ( fR) ] 


Figure  11a.  Typical  event  distributions  in  the  mt>(f)  plane 
for  event  data  that  is  not  corrected  for  noise 
contamination.  Noise  pulses,  when  viewed  in 
this  space  appear  roughly  as  shown  and  affect 
explosion  event  m^Cf^)  values  most  strongly, 
causing  population  overlap  at  low  magnitudes. 
The  population  boundaries  for  noise  and  events 
are  somewhat  source  and  receiver  dependent  due 
to  earth  structure  variations. 
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earthquake  limit  line 
(Parallel  to  discrimi¬ 
nation  line)  . 


NOISE  CORRECTED 
EARTHQUAKE 
POPULATION  BOUNDARY 


OBSERVED  MEAN 
VALUE  NOISE 
POPULATION 
BOUNDARY 


DISCRIMINATION 

LINE 


EXPLOSION 
LIMIT  LINE 
(Parallel 
to  discrimi¬ 
nation  line) 


NOISE  CORRECTED 
EXPLOSION 

POPULATION  BOUNDARY 


"CERTAIN"  EXPLOSION  EVENT 
IDENTIFICATION 


NOISE 


POINT  "UNCERTAINTY  ELLIPSE* 


NOISE  CORRECTED  EVENT  POINT,  WITH 
ERROR  BOUNDS  BASED  ON  UNCERTAINTY 
IN  THE  DETERMINISTIC  NOISE 
CORRECTION 


HIGH  FREQUENCY  BODY  WAVE  MAGNITUDE 


Figure  lib.  Typical  event  distributions  in  the  m^Cf)  plane  for 
noise  corrected  event  data.  A  discrimination  line 
can  be  defined  on  the  basis  of  the  definition  of 
these  populations  using  known  events,  or  on  the 
basis  of  theoretical  predictions. 
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III.  VFM  DISCRIMINATION  RESULTS  FROM 
A  TEN  STATION  NETWORK 


3.1  VFM  RESULTS  FROM  TELESEISMIC  STATIONS 

The  m^f)  results  obtained  from  the  ten  stations  in¬ 
cluded  in  this  phase  of  the  experiment  demonstrate  that  the 
VFM  approach  can  discriminate  between  Eurasian  earthquakes 
and  explosions  recorded  at  both  teleseismic  and  regional 
distances.  Individual  station  performance,  however,  is 
largely  controlled  by  the  attenuation  properties  of  the  upper 
mantle  beneath  a  particular  station.  For  instance,  a  station 
overlying  high  Q  upper  mantle  (as  evidenced  either  by  the  siting 
of  the  station  in  a  shield  region  and/or  a  region  characterized 
by  large  negative  P  wave  travel-time  residuals)  provides  much 
greater  separation  of  earthquakes  and  explosions  in  the  m^ff) 
plane  representation  than  a  station  located  over  a  low  Q  upper 
mantle.  As  examples,  we  emphasize  results  from  two  stations  for 
which  a  variety  of  regional  geophysical  information  pertaining 
to  the  structure  of  the  upper  mantle  is  available.  The  two 
stations  are  Red  Lake,  Ontario  (RKON)  and  Albuquerque,  New 
Mexico  (ANMO) . 

Figure  12  shows  discrimination  results  obtained  from 
the  vertical  component  seismograms  of  short-period  P  waves 
recorded  at  RKON.  The  epicentral  distance  range  of  all  the 
Eurasian  events  to  RKON  is  53  to  103  degrees.  The  particular 
pair  of  m.  (f)  frequencies,  0.6  Hz  and  3.25  Hz,  yielded  the 
greatest  separation  of  the  majority  of  earthquakes  and  explo¬ 
sions  recorded  at  this  location. 

As  seen  in  Figure  12,  the  VFM  technique  discriminates 
very  effectively  between  earthquakes  and  explosions  over  the 
entire  magnitude  range  of  these  events.  Out  of  a  total  of  50 
shallow  earthquakes  (open  circles) ,  only  one  event,  number  34, 
plots  in  the  explosion  population.  This  particular  earthquake 


47 


SYSTEMS  SCIENCE  AND  SOFTWARE 


o 


a  s 


O'  • 
vo  n 
*-»  w 


<? 


o 

o 


_ o  a5 
5< 


a  X-®5  ° 

B  o  r°o 

<5o  £CO  Z*3 
0-0 


Ow  00 


0 

O  •  OS 

o  o 

0 


XI 


B 

JC 

O  o 
•H  u"> 
H  , 
10  vl 
■U 

c  x: 

o  — 

-  tn 
<D  (U 
.*  .* 
10  (0 
i-5  3 
O’ 

•a  x: 

<V  44 
OS  >4 
(0 

4J  d) 

<n 

3 
O 


T3 

Q)  i-H 
4J  r-t 
(0  (0 
O  43 

O  to 


C  0 
O  -rt 

•h  a 
xj  at 
(0  T3 

4J 

to  to 

a> 

<U  r-4 

j=  y 

44  >4 
•H 
ti  y 
0 

44  C 

a> 

to  a 

.u  o 

rH 

3 

to  . 
0)  — 
U  2 
O 
S  * 
fc  a 
>  — 


a)  *— i 
mh  a 
aj  ct>  x 

M  C  V 
(0  <0 
3  -H  >0 

tr  vx  <u 
to  -u  g 
3 

c  c  to 
<u  a)  <u 
a  a  w 
o  o  a 


!=H  9  ’0) 


tU 

u 

3 

C> 


48 


SYSTEMS  SCIENCE  AND  SOFTWARE 


intermediate  focal  depths  (50  <  h  <  250  km), 
s  deep  events  (h  >  250  km)  and  closed  circles 


is  one  of  several  aftershocks  following  a  large,  m^  =  5.7, 
Tibetan  earthquake  (Event  Number  30).  The  remaining  three 
earthquakes  that  fail  to  discriminate  at  RKON  are  all  rela¬ 
tively  deep  focus  events.  The  focal  depths  reported  in  the 
Preliminary  Determination  of  Epicenters,  Monthly  Listings, 
published  by  the  U.  S.  Geological  Survey  (USGS)  are  as  follows: 
Event  28,  157  km;  Event  147,  479  km;  Event  169,  118  km.  As 
noted  in  previous  studies  (Savino  and  Archambeau,  1974;  Savino, 
et  al. ,  1975),  the  majority  of  events  that  fail  to  discriminate 
using  the  VFM  technique  are  deep.  We  point  out,  however,  that 
while  these  three  deep  earthquakes  plot  in  the  explosion  popu¬ 
lation  in  Figure  12,  16  other  deep  events,  including  two  with 
reported  depths  in  excess  of  500  km,  plot  well  within  the 
shallow  earthquake  population. 

Other  than  Event  number  79,  a  small  presumed  explosion 
in  eastern  Kazakhstan,  the  remaining  nine  explosions  recorded 
at  RKON  are  well  separated  from  the  earthquakes  in  Figure  12. 

It  is  particularly  interesting  to  note  that  while  there  are 
only  ten  explosions  recorded  at  this  station,  these  ten  events 
originate  in  four  separate  Eurasian  source  regions.  These 
include  Novaya  Zemlya,  Lake  Baikal,  eastern  Kazakhstan  and 
north  of  the  Caspian  Sea. 

As  for  the  attenuation  properties  of  the  upper  mantle 
beneath  RKON  we  note  that  this  station  is  located  on  the 
Canadian  Shield  in  a  region  characterized  by  relatively  low 
heat  flow  and  large  negative  travel- time  residuals  (Masso, 
et  al. ,  1978).  The  negative  travel-time  residuals  for  tele- 
seismic  events  imply  relatively  high  velocities  in  the  deep 
crust  and  upper  mantle  beneath  RKON.  These  geophysical  data 
sets,  together  with  results  from  several  other  studies  (e.g., 
Brune  and  Dorman,  1963) ,  suggest  that  the  upper  mantle  in  this 
region  is  relatively  high  Q  (low  attenuation) .  In  a  later 
section  of  this  report,  we  attempt  to  demonstrate  how  the 
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interplay  of  a  high  Q  upper  mantle  and  the  prevailing  spectrum 
of  background  earth  noise  determine  the  performance  of  the  VFM 
discriminant  at  a  particular  station. 

VFM  discrimination  results  similar  to  those  obtained 
for  RKON  were  found  for  two  other  stations:  Kabul,  Afghanistan 
(KAAO)  and  Charters  Towers,  Australia  (CTAO) .  The  results  for 
CTAO  are  shown  in  Figure  13. 

While  the  number  of  available  event  recordings  from  CTAO 
(Figure  13)  is  rather  small,  the  general  separation  of  shallow 
earthquakes  and  presumed  explosions  is  quite  striking,  more 
than  0.5  units  for  either  a  fixed  low  or  high  frequency 
m^Cf)  level.  The  only  earthquake  that  plots  in  with  the  pre¬ 
sumed  explosions  is  an  event  with  intermediate  focal  depth 
(Event  169) .  The  presumed  explosion.  Event  number  21,  plotting 
in  with  the  shallow  earthquake  population  is  anomalous  at 
several  other  stations  and  is  addressed  in  more  detail  in 
Savino,  et  al.  (1980) . 

In  contrast  to  an  RKON- type  station,  we  have  the  VFM 
discrimination  results  for  ANMO  (Figure  14) .  This  station  is 
located  in  the  Basin  and  Range  Province  of  the  western  United 
States  and  is  in  a  region  characterized  by  large  positive 
travel-time  residuals  and  high  heat  flow  (Masso,  et  al. , 

1978) .  Numerous  geophysical  studies  performed  in  this  region 
(Thompson  and  Burke,  1974)  indicate  the  existence  of  an  ex¬ 
tensive  low  Q-low  velocity  zone  in  the  upper  mantle.  In  terms 
of  discrimination,  an  ( f )  plot  obtained  for  the  Eurasian 
events  recorded  at  ANMO  can  best  be  described  as  a  shotgun 
pattern  over  most  of  the  event  magnitude  range,  with  signifi¬ 
cant  overlap  of  explosions  and  earthquakes  of  all  focal  depths. 

A  detailed  comparison  of  the  VFM  results  for  ANMO 
(Figure  14)  with  those  for  RKON  (Figure  12)  reveals  some  inter¬ 
esting  features.  For  instance,  note  in  Figure  14  that  the  high 
frequency  mb(f)  estimates  at  ANMO  do  not  exceed  a  value  of 
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Figure  13.  VFM  results  for  the  station  located  at  Chart 
Towers,  Australia  (CTAO) .  Legend  for  syrobol 
sane  as  in  Figure  12. 
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Figure  14.  VFM  results  for  the  station  located  at 
Albuquerque,  New  Mexico  (ANMO) .  Legend 
for  symbols  same  as  in  Figure  12. 
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about  4.8.  This  is  for  a  frequency  of  2.25  Hz.  At  RKON 
(Figure  12),  however,  the  high  frequency  mfa(f)  estimates,  now 
at  3.25  Hz,  attain  values  as  high  as  5.2.  In  more  detail,  con¬ 
sider  the  four  explosions  (Events  19,  20,  22  and  81)  which  are 
recorded  at  both  stations.  The  high  frequency  mfa(f)  estimates 
for  these  events  are  between  0.5  and  1.3  units  higher  at  RKON 
than  at  ANMO.  The  low  frequency  estimates,  on  the  other  hand, 
differ  by  only  0.2  to  0.3  magnitude  units.  The  fact  that  two 
of  these  explosions  (Events  20  and  81)  are  located  in  eastern 
Kazakhstan,  one  (Event  19)  at  Novaya  Zemlya  and  the  fourth 
(Event  22)  north  of  the  Caspian  Sea  argues  against  a  systematic 
source  region  effect  being  responsible  for  differences  in  the 
high  frequency  spectral  levels  observed  at  these  two  stations. 
In  addition,  the  overlap  in  the  ranges  of  epicentral  distances 
of  the  common  events  in  Figures  12  and  14,  54  to  79  degrees  at 
RKON  and  71  to  95  degrees  to  ANMO,  argues  against  possible 
variations  in  the  lower  mantle,  where  the  rays  bottom,  as  the 
cause  of  the  observed  spectral  differences.  Given  the  obser¬ 
vations,  the  most  likely  explanation  is  that  the  known  dif¬ 
ferences  in  high  frequency  attenuation  beneath  these  two 
recording  sites  are  responsible  for  the  relative  positions  of 
the  Eurasian  events  in  the  m^tf)  planes  shown  in  Figures  12 
and  14. 

Other  stations  which  behave  like  ANMO,  and  in  some  cases 
are  known  to  be  situated  in  regions  characterized  by  large 
positive  travel- time  residuals,  are  the  following:  Mashhad, 

Iran  (MAIO) ;  Chiang  Mai,  Thailand  (CHTO) ;  Taipei,  Taiwan  (TATO) ; 
and  Zongo  Valley,  Bolivia  (ZOBO) .  Figure  15  shows  VFM  results 
from  three  stations  in  the  ANMO  category.  These  stations  are 
TATO  (Figure  15a) ,  CHTO  (Figure  15b)  and  ZOBO  (Figure  15c) . 
While  the  number  of  events  available  for  analysis  at  TATO  and 
ZOBO  is  fairly  small,  in  general,  the  dominant  feature  of  these 
three  station  m^ ( f )  plots  is  the  lack  of  separation  of  the 
earthquake  and  explosion  populations. 

53 


SYSTEMS.  SCIENCE  A NO  SOFTWARE 


(0.4  HZ) 


TATO 


(2.0  Hz) 


'Figure  15a.  VFM  results  for  all  events  recorded  at  the 
station  in  Taipei,  Taiwan  (TATO) .  Legend 
for  symbols  same  as  in  Figure  12. 
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.gure  15c.  VFM  results  for  the  station  located  at  Zongo 
Valley,  Bolivia  (ZOBO) .  Legend  for  symbols 
same  as  in  Figure  12. 
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The  remaining  stations  at  teleseismic  distances  from 
the  Eurasian  events  are  Houlton,  Maine  (HNME)  and  the  center 
subarray  at  LASA  (LAO) .  HNME  is  located  in  the  Appalachian 
Highlands  and  LAO  in  the  Great  Plains  Geologic  Province. 

These  two  regions  are  characterized  by  essentially  average 
values  of  heat  flow,  Pn  velocities  and  crustal  thicknesses 
as  compared  to  either  the  Canadian  Shield  (RKON)  or  Basin  and 
Range  (ANMO)  Provinces.  As  seen  in  Figures  16a  and  16b,  the 
VFM  discrimination  results  for  HNME  and  LAO  are  intermediate 
to  those  for  RKON  and  ANMO. 


3.2  VFM  RESULTS  FROM  STATIONS  AT  REGIONAL  AND  TELESEISMIC 

DISTANCES 

In  the  following,  we  summarize  the  most  important  re¬ 
sults  obtained  from  three  stations  that  recorded  events  at 
both  regional  and  teleseismic  distances.  These  stations  are 
KAAO,  ILPA  and  MAIO. 

Figure  17  shows  the  locations  of  18  shallow  earthquakes 
(open  circles)  and  ten  presumed  explosions  (closed  circles) 
relative  to  the  Seismic  Research  Observatory  (SRO)  station, 
KAAO.  The  epicentral  distance  range  of  these  events  to  KAAO 
is  15  to  23  degrees,  with  the  nine  presumed  explosions  in 
eastern  Kazakhstan  at  approximately  17  degrees.  The  VFM  re¬ 
sults  for  these  28  events  are  given  in  Figure  18a.  As  can  be 
seen,  the  presumed  explosions  clearly  separate  from  the  shallow 
earthquakes  over  the  entire  magnitude  range.  In  particular, 
the  VFM  technique  successfully  discriminates  Events  80,  162, 

163  and  272,  all  of  which  are  close  to  explosion  test  sites. 

The  implication  here  is  that  differences  in  the  propagation 
path  are  minimized  and  we  are  seeing  differences  in  explosion 
and  earthquake  source  spectra. 

A  very  important  point  to  note  about  the  m^(f )  estimates 
in  Figure  18a  is  the  value  of  the  filter  center  frequency  used 
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VFM  results  for  the  station  located  at 
Houlton,  Maine  (HNME) .  Legend  for  sym¬ 
bols  same  as  in  Figure  12.  We  suspect 
that  the  digital  gain  factors  supplied 
by  the  SDAC  for  the  eight  shallow  earth¬ 
quakes  with  mb  (0.5  Hz)  values  of  6.5 
and  greater  are  incorrect  -  too  high, 
in  fact.  Incorrect  gain  factors,  how¬ 
ever,  only  move  event  points  along  diag¬ 
onals  in  the  n^if)  plane  and  do  not 
impair  separation  of  event  populations. 
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Figure  16b.  VFM  results  for  the  center  subarray  at  LASA 
(LAO) .  Legend  for  symbols  same  as  in  Figure 
12. 
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Figure  17.  Map  of  Eurasian  earthquakes  (open  circles)  and 
presumed  explosions  (closed  circles)  located 
between  15  and  23  degrees  from  the  SRO  station 
at  Kabul,  Afghanistan  (KAAO) . 
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Figure  18a.  VFM  discrimination  results  for  Kabul, 
Afghanistan  (KAAO)  for  the  earthquakes 
(open  circles)  and  presumed  explosions 
(closed  circles)  mapped  in  Figure  17. 

These  events  are  located  between  15° 
and  23°  from  KAAO.  Note  the  separation 
between  Event  272  and  the  nearby  nine 
presumed  explosions  at  eastern  Kazakhstan. 


! 
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for  the  high  frequency  estimates  (i.e.,  4.0  Hz).  This  is  quite 
high  compared  to  previous  results  (Savino,  et  al. ,  1979)  and  is 
probably  a  result  of  propagation  over  a  fairly  high  Q  path  in 
the  topmost  mantle  considering  the  far  regional  distance  range 
of  the  events.  The  results  in  Figure  18b  support  this  inter¬ 
pretation.  In  this  figure,  we  give  m^tf)  results  for  events 
at  teleseismic  distances  from  RAAO.  The  particular  pair  of 
high  and  low  frequencies  (2.75  and  0.4  Hz,  respectively)  in 
this  figure  gave  the  best  discrimination  results.  Most  of  the 
events  in  Figure  18b  that  fail  to  discriminate  can  be  explained 
in  terms  of  either  large  focal  depth  (the  open  triangles  and 
squares)  or  anomalous  attenuation  in  the  upper  mantle  near  the 
source  region  of  two  of  the  presumed  explosions,  21  and  16 
(Savino,  et  al . ,  1980).  The  point  we  want  to  note  here  is  that 
the  frequency  corresponding  to  the  high  frequency  m,  (f)  esti¬ 
mate  is  lower  in  the  case  of  events  at  teleseismic  distances 
than  at  far  regional  distances. 

A  similar  dependence  of  the  frequency  value  of  the  high 
frequency  m^ff)  estimates  on  event  distances  is  shown  in  Figures 
19a  and  19b.  The  events  in  question  occurred  throughout  the 
same  general  region  shown  in  the  map  in  Figure  17.  The  record¬ 
ing  station  in  this  case  is  the  Iranian  Long  Period  Array  (ILPA) . 
In  Figure  19a  we  see  an  interesting  separation  of  (1)  earth¬ 
quakes  and  the  only  presumed  explosion  recorded  in  this  distance 
range  at  ILPA,  and  (2)  earthquakes  at  distances  of  6.9  to  9.2 
degrees  from  earthquakes  at  distances  between  .2.5  and  19 
degrees.  The  presumed  explosion.  Event  number  22,  is  approxi¬ 
mately  12.5  degrees  from  ILPA.  The  VFM  results  in  Figure  19a 
are  primarily  based  on  Pg  for  the  regional  events  included 
within  the  solid  line  and  mantle  P  waves  for  all  the  other 
events.  While  the  presumed  explosion  separates  extremely  well 
from  events  at  equal  or  greater  distances,  the  separation 
between  it  and  the  six  earthquakes  at  closer  distances  is  much 
less.  This  is  most  likely  due  to  the  more  efficient  high  fre¬ 
quency  Pg  propagation  in  the  case  of  the  earthquakes. 
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Figure  19a.  VFM  results  for  regional  events  recorded  at 
ILPA.  The  dashed  and  solid  curves  group 
events  in  the  two  distance  ranges.  The  epi- 
central  distance  of  the  presumed  explosion. 
Event  22,  to  ILPA  is  12.5  degrees.  See 
Figures  12  and  18b  for  definitions  of  symbols. 
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Figure  19b.  VFM  results  for  teleseismic  events  recorded 
at  ILPA.  See  Figures  12  and  18b  for  defini¬ 
tions  of  symbols. 
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Figure  19b  shows  VFM  discrimination  results  for  events 
at  telseismic  distances  (A  >_  30°)  from  ILPA.  As  we  saw  for 
KAAO,  going  from  regional  to  teleseismic  distances  results  in 
a  reduction  in  the  frequency  of  the  high  frequency  mfa(f)  esti¬ 
mates  that  provide  optimal  discrimination. 

In  Figure  20a  we  give  the  VFM  results  for  the  SRO  sta¬ 
tion  located  in  Mashhad,  Iran  (MAIO)  for  most  of  the  events 
mapped  in  Figure  17.  The  first  point  to  be  noted  is  that 
while  the  distance  range  to  MAIO  (i.e.,  18.7°  <_  A  <  28.7°)  is 
comparable  to  KAAO  (15°  <_  A  £  23°,  Figure  18a),  the  highest 
frequencies  reliably  observed  at  MAIO  are  only  2.25  Hz. 
Secondly,  it  is  immediately  obvious  that  the  degree  of  separa¬ 
tion  of  event  populations  at  MAIO  is  significantly  less  than 
the  separation  seen  at  KAAO. 


For  events  at  teleseismic  distances  from  MAIO  we  have 
the  results  plotted  in  Figure  20b.  The  epicentral  distances 
of  all  the  events  are  greater  than  30  degrees.  As  mentioned 
previously,  MAIO  behaves  like  ANMO  in  terms  of  VFM  type  dis¬ 
crimination,  yielding  essentially  no  separation  of  earthquake 
and  explosion  populations. 

Differences  in  the  performance  of  the  VFM  discrimination 
technique  at  KAAO,  ILPA  and  MAIO  are  primarily  attributed  to 
differences  in  attenuation  in  the  upper  mantle  beneath  these 
stations.  Support  for  this  statement  comes  from  a  study 
(Masso,  et  al . ,  1978)  of  teleseismic  P-wave  residuals  deter¬ 
mined  for  a  global  network  of  seismograph  stations,  including 
several  stations  located  within  the  region  being  considered 
here.  Figure  21  is  a  map  of  P-wave  residuals  determined  from 
bulletins  of  the  International  Seismological  Center  (ISC)  for 
stations  (solid  triangles)  in  this  region.  Each  residual  value 
plotted  in  this  figure  is  based  on  a  minimum  of  100  globally 
distributed  events.  What  we  observe  is  that  stations  in  the 
vicinity  of  KAAO,  including  the  KAAO  location  itself,  are 
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Figure  20a.  VFM  results  for  the  station  at  Mashhad, 
Iran  (MAIO)  for  events  in  the  distance 
range  18.7°  <  &  £  28.7°.  The  open 
circles  are  earthquakes,  the  closed 
circles  presumed  explosions. 
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Teleseismic  P-wave  residuals  based  on  data 
from  the  bulletins  of  the  ISC  for  seismo¬ 
graph  stations  (closed  triangles)  in  the 
region  of,  and  including,  MAIO,  KAAO  and 
ILPA.  These  latter  three  stations  are 
denoted  by  the  closed  circles. 
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characterized  by  large  negative  residuals.  MAIO  and  nearby 
stations,  on  the  other  hand,  exhibit  large  positive  residuals. 
While  we  lack  information  from  ILPA,  we  note  that  this  station 
is  surrounded  by  stations  to  the  north  with  relatively  large 
positive  residuals  and  one  station  to  the  southwest  with  a 
negative  residual.  Based  on  similar  measurements  for  stations 
in  other  parts  of  the  world,  a  regional  distribution  of  large 
negative  P-wave  residuals  (i.e.,  near  KAAO)  is  attributed  to 
a  relatively  fast,  low-attenuating  upper  mantle  whereas  large 
positive  residuals  (e.g.,  near  MAIO)  are  usually  associated 
with  a  slow,  highly  attenuating  upper  mantle. 
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IV.  DISCUSSION  AND  CONCLUSIONS 


In  the  previous  sections  of  this  report  we  introduced 
arguments  to  explain  various  characteristics  of  the  VFM  dis¬ 
criminant  observed  at  the  different  stations  included  in  this 
experiment.  In  particular,  we  suggested  that  the  combined 
effects  of  background  earth  noise  and  anelastic  attenuation 
play  an  important  role  in  determining  the  performance  of  the 
VFM  technique  at  the  individual  stations.  In  the  following 
we  will  develop  these  arguments  in  more  detail. 


4.1  INFLUENCE  OF  PATH  AND  SOURCE  EFFECTS  ON  VFM  DISCRIMINATION 

The  way  in  which  we  view  anelastic  attenuation  in  the 
upper  mantle  and  background  earth  noise  at  a  particular  record¬ 
ing  location  contributing  to  VFM  discrimination  is  as  follows. 
Anelastic  attenuation  beneath  a  given  station  affects  P  waves 
from  both  earthquakes  and  explosions  at  teleseismic  distances 
in  the  same  way,  basically  imposing  an  exponentially  shaped 
low  pass  filter  on  the  incoming  signals.  Both  event  popula¬ 
tions,  when  viewed  in  an  m^(f)  plane,  are  supposedly  shifted 
by  the  same  amount  with  the  largest  changes  being  in  the  high 
frequency  m^ff)  values  (i.e.,  recall  the  behavior  of  explosions 
recorded  at  RKON  relative  to  ANMO ,  in  Figures  12  and  14,  re¬ 
spectively,  in  Section  III) .  However,  because  of  the  presence 
of  a  station  noise  background  that  is  largely  unaffected  by 
the  upper  mantle,  the  event  populations  cannot  freely  shift 
through  the  noise  population  in  the  m^ff)  plane  (Savino,  et  al . , 
1979).  Thus,  earthquakes  (which  occupy  a  region  in  the  m^ff) 
plane  close  to,  and  even  overlapping  in  the  case  of  small  mag¬ 
nitude  events,  the  noise  population)  will  approach  and  eventually 
totally  merge  with  the  noise  population  as  the  station  upper 
mantle  Q  becomes  lower  and  lower.  When  the  earthquake  m^ff) 
estimates  are  within  or  very  close  to  the  noise  population  they 
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can  no  longer  move  in  the  m^f)  plane  and  so  will  not  shift  any 
more  even  if  the  station  Q  is  further  reduced. 

The  explosion  population,  as  a  function  of  lower  and 
lower  mantle  Q,  will  also  appear  to  move  toward  the  noise 
population  in  the  same  way  as  the  earthquakes,  but  because 
they  occupy  a  region  with  larger  high  frequency  m^f)  values 
than  those  for  earthquakes  (for  whatever  reason) ,  they  will 
continue  to  shift  toward  the  earthquake  population  with  lower¬ 
ing  station  Q,  even  then  the  earthquakes  are  no  longer  able  to 
shift.  In  this  case,  when  the  earthquake  population  reaches 
the  "noise  stop"  then  the  explosions  will  begin  to  converge  on 
the  earthquakes  with  a  resultant  reduction  in  the  "discrimina¬ 
tion  gap"  seen  at  high  Q  stations.  Thus,  for  stations  in  very 
low  mantle  Q  zones  and/or  with  high  noise  levels,  one  would 
expect  to  see  converging  or  overlapping  event  populations  in 
the  mb(f)  plane.  This  appears,  in  fact,  to  be  the  case  when 
m^Cf)  discrimination  at  high  Q  mantle  stations  ( RKON ,  KAAO) 
is  compared  to  mfa(f)  discrimination  at  stations  over  low  Q 
mantle  zones  (ANMO,  MAIO) .  Also,  especially  for  event  popula¬ 
tions  defined  by  small  numbers  of  events,  this  would  cause 
variations  in  the  position  and  slope  of  the  discrimination 
lines  separating  earthquakes  and  explosions  at  stations  over 
different  mantle  Q  structures.  Since  high  frequencies  are 
used,  this  variation  is  quite  large. 

The  sensitivity  of  the  VFM  technique  to  variations  in 
the  propagation  path  raises  another  possibility  concerning  the 
reasons  (source  versus  path  effects)  for  its  overall  applica¬ 
bility.  For  instance,  it  might  be  argued  that  the  compressional 
wave  radiation  spectra  from  earthquakes  and  explosions  are  not 
different  in  any  significant  way.  Rather,  the  reason  for  the 
discrimination  observed  in  Section  III  is  that  the  Eurasian 
explosions  are  shallow  events,  so  there  is  low  frequency  can¬ 
cellation  due  to  pP  interference.  In  addition,  since  they 
occur  in  regions  that  are  aseismic  and,  in  several  cases 
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tectonically  stable  with  shield-like  upper  mantle  structures, 
the  paths  in  the  source  region  are  "high  Q  paths."  On  the 
other  hand,  all  the  earthquakes  are  reasonably  deep  with  little 
or  no  pP  low  frequency  cancellation  and  occur  in  tectonically 
active  zones  within  which  "low  Q  mantle"  paths  can  be  expected. 
Thus,  the  differences  in  the  observed  spectral  content,  as 
manifested  by  the  m^ ( f )  discrimination  data  for  the  Eurasian 
events,  are  due  to  differences  in  both  low  frequency  pP  can¬ 
cellation  and  attenuation  between  events  in  different  tectonic 
regions.  In  this  case,  the  high  frequencies  from  the  earth¬ 
quakes  would  be  expected  to  be  more  strongly  reduced  by  attenua¬ 
tion  due  to  their  occurrence  in  low  Q  mantle  regions.  While 
this  argument  could  explain  the  overall  separation  of  earth¬ 
quakes  and  explosions  obtained  at  the  high  Q  stations  (e.g., 
RKON,  KAAO,  CTAO) ,  when  we  consider  the  results  in  more  detail 
we  find  that  this  argument  encounters  several  difficulties. 

In  order  to  address  these  difficulties  we  repeat  the  VFM 
discrimination  results  for  KAAO  and  RKON  in  Figure  22.  Two 
tests  that  we  can  perform,  given  the  discrimination  results 
from  this  experiments,  for  distinguishing  between  source  dif¬ 
ferences  or  propagation  path  (Q)  differences  as  the  underlying 
reason  for  m^Cf)  discrimination  are  (1)  to  compare  m^Cf)  re- 
suits  for  collocated  earthquakes  and  explosions  and  (2)  exa¬ 
mine  the  mb(f)  behavior  of  earthquakes  for  which  the  entire 
propagation  paths  are  presumably  high  Q  —  certain  shallow 
trench  events  and  deep  earthquakes  recorded  at  Class  I  stations. 

With  respect  to  the  first  test,  we  refer  the  reader  to 
Figure  17.  In  that  figure  we  mapped  the  locations  of  18  shallow 
earthquakes  in  the  regional  distance  range  from  KAAO.  Of  par¬ 
ticular  interest  is  Event  272,  a  shallow  earthquake,  located 
within  150  km  of  the  eastern  Kazakhstan  test  site.  As  seen  in 
Figure  22a,  m^ff)  estimates  for  Event  272  plot  within  the 
population  of  earthquakes  recorded  at  KAAO,  with  a  difference 
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Figure  22a.  VFM  discrimination  results  for  Kabul, 
Afghanistan  (KAAO)  for  the  earthquakes 
(open  circles)  and  presumed  explosions 
(closed  circles)  mapped  in  Figure  17. 

These  events  are  located  between  15° 
and  23°  from  KAAO.  Note  the  separation 
between  Event  272  and  the  nearby  nine 
presumed  explosions  at  eastern  Kazakhstan. 
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(3.25  Hz) 


Figure  22b.  VFM  results  for  the  station  located  at  Red  Lake, 
Ontario  (RKON) .  This  figure  is  similar  to 
Figure  12  except  that  several  additional  earth¬ 
quakes  in  the  Japan-Kuril  arcs  and  Hindu  Kush 
regions  have  been  identified. 
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of  at  least  one  m^  unit  in  the  high  frequency  magnitude  esti¬ 
mates  for  this  event  compared  to  those  for  presumed  explosions 
at  eastern  Kazakhstan.  Similar  comments  apply  to  the  three 
shallow  earthquakes  west  of  the  Caspian  Sea  (Events  80,  162, 

163)  compared  to  the  presumed  explosion  in  that  region  (Event 
22) ,  although  these  events  are  not  nearly  as  close  as  Event 
272  is  to  the  nine  presumed  explosions  at  eastern  Kazakhstan. 

The  conclusion  we  draw  from  this  comparison  is  that,  with  dif¬ 
ferences  in  propagation  paths  minimized,  the  reason  the  earth¬ 
quakes  and  presumed  explosions  separate  in  the  m^tf)  plane 
(Figure  22a)  is  the  enriched  high  frequency  source  spectra  of 
the  presumed  explosions  relative  to  the  earthquakes. 

For  the  second  test,  we  select  earthquakes  and  recording 
stations  for  which  the  entire  paths  to  the  stations  are  most 
likely  high  Q  paths.  Event-station  combinations  that  most 
likely  satisfy  this  criterion  are  (1)  shallow  trench  earthquakes 
with  receiving  stations  on  the  oceanic  side  of  the  events  and 
(2)  deep  earthquakes  (i.e.,  below  low  Q  zones)  recorded  at 
stations  sited  on  shields.  Given  these  experimental  conditions, 
if  the  earthquake  and  explosion  source  spectra  are  really  not 
much  different,  then  all  the  explosions  should  appear  in  the 
earthquake  populations. 

The  m^tf)  results  from  all  the  Class  I  stations  for  events 
along  the  Japan-Kuril-Kamchatka  arcs  and  in  the  Hindu  Kush 
region  argue  strongly  for  differences  in  source  spectra  —  see 
Tables  A. 1  and  A. 2  in  Appendix  A  of  the  report  by  Savino,  et  al. 
(1980)  for  Events  47,  48,  49,  50,  62,  63,  64,  69,  143,  144,  145, 
147,  165  and  169.  As  an  example,  we  repeat  m^ff)  results  for 
the  Class  I  station  RKON  in  Figure  22b.  Several  additional 
shallow  and  deep  events  occurring  along  the  island  arc  systems 
and  in  the  Hindu  Kush  region  are  identified  on  this  figure. 

First  we  note  that  the  only  shallow  earthquake  that  appears 
explosion-like  at  this  station  is  Event  34,  an  aftershock  of 
a  large  Tibetan  earthquake  (Event  30) .  None  of  the  shallow 
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trench  earthquakes  (e.g.,  Events  47  and  50)  plots  in  the  explo¬ 
sion  population.  The  second  point  is  that  out  of  a  total  of 
19  deep  events  analyzed  at  this  station,  only  three  appear 
explosion-like.  The  fact  that  the  locations  of  these  three 
events  (147  and  169  in  the  Kuril  Islands,  28  in  the  Hindu  Kush) 
are  common  to  the  16  events  that  appear  earthquake- like ,  argues 
against  systematic  path  effects  being  responsible  for  discrimina¬ 
tion.  Rather,  based  on  the  multistation  results  for  all  these 
events  (Savino,  et  al. ,  1980) ,  the  preferred  explanation  is  that 
the  source  spectra  for  the  earthquakes  and  explosions  are  dif¬ 
ferent  and  the  occasional  explosion-like  behavior  of  certain 
earthquakes  at  individual  stations  can  be  attributed  to  azimuthal 
variations  in  source  radiation  patterns. 


4.2  CONCLUSIONS 

The  principal  conclusions  resulting  from  the  discrimina¬ 
tion  experiment  are  the  following: 

•  VFM  discrimination  can  be  achieved  to  m.  'v  4 
at  Class  I  regional  stations. 

•  VFM  discrim  .nation  based  on  both  regional  and 
teleseismic  stations  is  approximately  95  per¬ 
cent  effective  for  all  the  events  in  this 
experiment. 

•  Stations  can  be  classified  as  to  VFM  dis¬ 
crimination  effectiveness/reliability  using 
P-wave  travel-time  residual  data  and  other 
information  pertaining  to  the  attenuation 
properties  of  the  upper  mantle  beneath  a 
site. 

•The  VFM  technique  is  an  effective  discriminant 
principally  because  of  the  enriched  high  fre¬ 
quency  source  spectra  of  explosions  relative 
to  earthquakes  with  the  same  spectral  level 
near  two  seconds  period. 
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APPENDIX  A 

TABULATION  OF  DISCRIMINATION  PARAMETERS 
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This  appendix  presents  a  complete  summary  of  the  fixed- 
frequency  discrimination  data  base  for  the  ten  stations  in¬ 
cluded  in  this  experiment:  CHTO,  CTAO,  KAAO,  MAIO,  TATO,  ZOBO, 
LAO,  ILPA,  RKON  and  HNME.  Included  in  the  tabulation  for  each 
station  are  only  those  individual  events  which  satisfy  the  dis¬ 
tance  range  criteria  defined  below  after  excluding  certain 
seismograms  judged  to  be  erroneous.  The  individual  seismograms 
were  first  processed  by  MARS,  with  the  primary  signal  detection 
and  noise  information  generated  for  each  of  40  narrow-band 
filters  being  saved.  All  of  these  output  for  approximately  513 
seismograms  were  merged  to  create  a  common  Multi-Frequency 
Discrimination  Data  Base  of  over  2.9  megabytes.  A  special  pur¬ 
pose  S3  Seismic  Data  Management  System  (SDMS)  was  designed  to 
operate  as  a  MARS  post-processor.  The  SDMS  was  applied  to 
this  common  data  base  for  all  ten  stations  first  to  select  all 
of  the  events  processed  for  a  given  station  and  to  generate  a 
new  single-station  multi-frequency  data  base  of  more  modest 
size?  a  set  of  single-station  data  bases  was  thereby  created 
for  more  efficient  post-processing.  Subsequent  applications 
of  the  SDMS  were  then  performed  in  order  (1)  to  exclude 
erroneous  seismogram  determined  by  prior  runs?  (2)  to  apply 
minimum  S/N  and  distance  range  criteria?  (3)  to  apply  the  VFM 
discrimination  algorithm  and  finally  (4)  to  select  an  appro¬ 
priate  set  of  discrimination  parameters  suitable  for  each 
individual  station.  Once  these  station-dependent  discrimina¬ 
tion  parameters  were  determined,  the  SDMS  was  used  to  create 
a  series  of  output  files  containing  the  final  discrimination 
results  for  each  station. 

These  output  files  for  the  set  of  stations  listed  above 
were  then  merged  to  c  eate  the  Fixed-Frequency  Discrimination 
Data  Base  for  the  stations  listed  above.  The  station  tabula¬ 
tions  given  here  summarize  these  results.  The  tabulation 
column  headings  are  described  below. 
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EVENT 


-  Event  number. 


DISTANCE  -  Epicentral  distance  in  degrees. 

FREQ  -  The  fixed  high  and  low  discrimination 
frequencies  fH  and  f^  (Hz)  . 


The  magnitude  estimate  at  the 
discrimination  frequencies  given 
by  the  relation, 

MB  (F )  =  log1Q  [A(F5  •  F]  +  b, 

where  A(F)  is  the  fitted  value  of 
the  noise  corrected  amplitude  at 
the  discrimination  frequency, 

?  =  f  or  f_ ,  and  b  is  distance 

ri  Li 

correction. 

MB (  +  )  ,  MB ( - )  -  Uncertainty  in  the  magnitude  esti¬ 
mate  as  given  by  the  relation, 

MB  (  +  )  =  log1Q  [  (A  (F)  ^iA^F))  •  F] 
+  b  -  MB  (F)  . 

Noise  Information 

AMP  -  The  fitted  value  for  the  noise  ampli¬ 

tude  at  the  fixed  high  (f^)  or  low 
(f_)  frequency.  These  fitted  values 
are  computed  from  the  least  squares 
polynomial,  of  degree  two  which  best 
approximates  the  unweighted  set  of 
log  A^Cf^)  values  over  the  frequency 
set  {fT>  or  {fH>,  where  A^(f^)  is  the 
mean  noise  amplitude  at  frequency  f.^. 
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MB  -  The  magnitude  estimate  corresponding  to 

the  fitted  value  of  the  noise  amplitude 
at  the  fixed  high  or  low  frequency. 


(S/N) **2  -  The  simple  mean  signal-to-noise  power 

_  2 

[A*  (f^) /A^  (f,^)  ]  averaged  over  the  high 
{ f H }  or  low  {fL>  frequency  bands. 

Signal  Arrival  Information 

TG  -  The  noise  weighted  mean  signal  group 

arrival  time  (seconds)  estimated  for 
the  high  or  low  frequency  bands  mea¬ 
sured  with  respect  to  the  event  origin 
time; 


t  -  S',  w. 
^  k-1 


k  yv 


/ 

£  w,  , 
k»l 


where  t^(f^)  is  the  identified  signal 
group  arrival  time  for  frequency  f,  , 
w,  is  the  weight  factor,  and  the  sum- 
mation  is  over  the  high  or  low  fre¬ 
quency  bands. 


TG3AR  -  Mean  group  arrival  time  t  measured 
from  the  event  origin. 


TSIGMA  -  Mean  deviation  in  t  ;  t  (first  line) , 
_  9  ~ 

t„  (second  line) . 

3 


NOISE 

WINDOW 


Length  in  seconds  of  the  noise  window 
used  in  the  calculation  of  mean  noise 
amplitudes. 


NFFT  -  The  number  of  frequency  samples  used 
in  the  discrete  Fourier  transform. 
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